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ABSTRACT 

In order to investigate the growth of super- massive black holes (SMBHs), we construct the black hole 
mass function (BHMF) and Eddington ratio distribution function (ERDF) of X-ray-sclcctcd broad- 
line AGNs at z ~ 1.4 in the Subaru XMM-Newton Deep Survey (SXDS) field. In this redshift range, a 
significant part of the accretion growth of SMBHs is thought to be taking place. Black hole masses of 
X-ray-selected broad-line AGNs are estimated using the width of the broad Mg II line and the 3000A 
monochromatic luminosity. Wc supplement the Mg II FWHM values with the Ha FWHM obtained 
from our NIR spectroscopic survey. Using the black hole masses of broad-line AGNs at redshifts 
between 1.18 and 1.68, the binned broad-line AGN BHMF and ERDF are calculated using the V mSiX 
method. To properly account for selection effects that impact the binned estimates, we derive the 
corrected broad-line AGN BHMF and ERDF by applying the Maximum Likelihood method, assuming 
that the ERDF is constant regardless of the black hole mass. We do not correct for the non-negligible 
uncertainties in virial BH mass estimates. If wc compare the corrected broad-line AGN BHMF with 
that in the local Universe, the corrected BHMF at z = 1.4 has a higher number density above 1O 8 M 
but a lower number density below that mass range. The evolution may be indicative of a down-sizing 
trend of accretion activity among the SMBH population. The evolution of broad-line AGN ERDF 
from z = 1.4 to indicates that the fraction of broad-line AGNs with accretion rate close to the 
Eddington-limit is higher at higher redshifts. 

Subject headings: galaxies: active — galaxies: statistics — galaxies: evolution — quasars: emission 
lines — quasars: general 



1. INTRODUCTION 

Since the discovery that super-massive black holes 
(SMBHs) sit at the centers of most massive galaxy in the 
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local Universe (e.g., IKormendv fc Richsto nc 1995|), de- 
termining their formation history remains one of big chal- 
lenges in astrophysics. It has been further determined 
that the mass of SMBHs correlates tightly with the phys- 
ical properties of their host spheroids (Mbh vs. Mb u i ge 
relations: e.g.. iMagorrian et al.l 119981 iMarconi fc Huntl 
[2001 IQultekin et al.N2009h . Such a correlation implies 
a physical connection between the growth histories of 
SMBHs and the sp heroidal components of galaxies (e.g. 
IBovle fc Terleviclj[l998h . 

Bolometric luminosities of AGNs reflect the mass ac- 
cretion rates of their SMBHs, therefore the luminos- 
ity function of AGNs and its cosmological evolution rc- 
flects the gro wth history of SMBHs through accretion 
(|Soltanlll982h . Cosmological evolution of AGN luminos- 



samples (e.g . lUeda et al.l 12003: ISilyerman et al 
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2008 




2011 



Croom et al . . 

Simpson et al.ll2012|) . The number density evolution of 
AGNs in different luminosity bins shows that higher 
luminosity AGNs, i.e. QSOs, have a peak at higher 
redshifts, the so-called "down-sizing" trend of cosmo- 
logical evolution of AGNs. The total amount of ac- 
creted matter estimated by integrating the luminosity 
functions over luminosity and redshift roughly matches 
to the estim ated mass density o f SMBHs in the lo- 
cal Universe (lYu fc Tremaind[200l IMarconi et al.ll2004t 
IShankar et al.H2009h . thus accretion is thought to be the 
dominant mode of SMBH growth. Appl ying the conti- 
nuity equation for the SMBH population, IMarconi et al.l 
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(12004ft evaluate the average growth curves of massive and 
less-massive SMBHs as a function of redshift. These re- 
sults imply that SMBHs grow rapidly at redshifts be- 
tween 1 and 2, and more massive SMBHs grow more 
than less massive SMBHs in the earlier Universe , as 
expected from the "down-sizing" trend of the AGN lu- 
minosity function. 

The luminosity of an AGN does not simply reflect the 
mass of its SMBH. In the calculations of the black hole 
growth history, the Eddington ratio, i.e. the ratio be- 
tween the observed accretion rate and the Eddington- 
limitcd accretion rate (A E dd), is assumed to be constant 
in AGNs with different luminosities and redshifts. As- 
suming a constant Eddington ratio for AGNs, their lumi- 
nosity directly corresponds to the mass of their SMBHs, 
and the mass-dependent growth history of SMBHs can 
be calculated. However, recent evaluations of the 
Eddington-ratio distribution of AGNs in the local uni- 
verse show that AGNs have a wide range of Eddington- 
ratio with no preferred value feauffma nn fc Heckmanl 
l2009t iSchulze fe Wisotzkl[20Tot) . Therefore, in order to 
quantitatively understand the accretion growth history 
of SMBHs, it is necessary to evaluate SMBH masses in 
any AGN sample for which the cosmological evolution of 
the luminosity function is evaluated. 

Tools are now available to measure black hole masses 
of broad-line AGNs. R everberation mapping of lo- 
cal broad-line AGNs (e.g.. IPeterson et al~ll2004l ) reveals 
the scaling relationship between the 5100A monochro- 
matic luminosity (L5100) of the broad -line AGN and the 
size o f its broad H/3 emitting region (|Kaspi et all I20001 
l2005f ). Utilizing this relationship, black hole masses of 
a large sample of broad-line AGNs can be estimated 
from their luminosities and broad H/ 3 line widths (Auh^) 
(e.g., IVestergaard k Peterson! [20061 ) from the relation- 
ship Mbh = /A w h^5ioo under the assumption that 
the b road-line region is virialized (Peterson k Wandell 
Il999t) . The factor / depends on the dynamical struc- 
ture of the broad-line region, and it is empirically de- 
termined by assuming that the black hole mass ob- 
tained from the reverberation mapping method and 
the velocity dispersion of its host bulge follow the 
Mb h vs. Ohnig p rela tion of local non-active galax- 
ies (jOnken et all |2004| ). This method has been ex- 
tended to black hole mass estimations using other broad 
lines such as Mg II AA2796,280 3A (jMcLure k Jarvis! 
[20021 IVestergaard fc Osmerll2009D. C IV AA 1548. 155lA 
( Vcstcrgaard 12002); IVestergaard k Peterson! 12009 ) . and 
Ha (|Greene fc Hdl2005D ."and using luminosities at other 
wavelen gths such as the 3000 A monochromatic lumi- 
nosity^ McLure fc Jarvisl I2002T ) , the Ha line luminos- 
ity llGreene k Holl2005D . and the hard X-ray luminosity 
l|Greene et al.ll2010al ). These relationships are calibrated 
against the black hole mass estimated by the reverber- 
ation mapping or that from the single-epoch broad-line 
width of H/3 and the monochromatic luminosity. Using 
these extended methods, black hole masses of broad-line 
AGNs at various redshifts can be estimated from their 
single-epoch optical spectra. 

Applying black hole mass estimates from single- 
epoch spectra to statistical samples of broad-line 
AGNs, black hole mass functions ( BHMFs) of broad- 
line AGNs can be evaluated (jWang et all 120061 : 
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Universe. ISchulze fc Wisotzkil (120101 ) derived the BHMF 



and Eddington-ratio distribution function (ERDF) of 
broad-line AGNs detected in the Hamburg/ESO AGN 
survey. They corrected the effects of the flux limits of 
their survey in their evaluation of the broad-line AGN 
BHMF and ERDF, i.e. the fact that the low-mass end 
of the sample only covers high Eddington ratio AGNs, 
by assuming that the ERDF is constant regardless of 
the black hole mass and applying Maximum likelihood 
method. Hereafter we label BHMF and ERDF derived 
by Vmax method with binned and those corrected for the 
detection limit by Maximum likelihood method with cor- 
rected. The corrected broad-line AGN BHMF covering 
Mbh values between 1O 6O M and 10 95 M Q and A E dd 
down to 0.01 shows rather steep decrease in number den- 
sity as a function of mass with no significant break in 
the mass range covered. Their corrected broad-line AGN 
ERDF shows a steep decline at the Eddington limit and 
a steep increase in the number density down to A E dd of 
0.01, following power-law with index of ~ —1.9. 

Cosmological evolution of the BHMFs of broad-line 
AGNs has also been examined using large samples of 
broad-line AGNs from the Sl oan Digital Sky Survey 
(SPSS) using the V m ^ meth od (jVestergaard et all [20081 : 
IVestergaard fc Osmerl 120091) and a Bayesian approach 
(IKellv et all 120101: IShen fc Kellvl 12012 IKellv fc Shen I 
20121). IKellv et all l20ioi) . IShen fc Kellvl (I2012D . and 
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of the BHMF of broad-line AGNs in the redshift range 
between 0.3 and 5 by applying a Bayesian approach 
(IKellv et al.ll2009l ). Hereafter we label BHMF and ERDF 
derived with the Bayesian approach with estimated. The 
estimated BHMF of broad-line AGNs shows an increase 
in number density above Mbh = 1 x 10 9 M q from z = 
to 2. In contrast, lower mass SMBHs show a relatively 
flat number density evolution up to z = 2. The "down- 
sizing" trend expected from the AGN luminosity function 
is confirmed by the steeper decrease of active SMBHs in 
the higher mass range from z = 2 to 0. However, it needs 
to be noted that broad-line AGNs in the SDSS sample 
only cover large Mbh and large A E dd AGNs in the red- 
shift range. For example, at z = 1.5 the sample is only 
30% complete dow n to Mbh ~ lO 9 Af and A E dd ~ 0.6 
(|Kellv et al.H2010T ). Due to the shallow detection limit, 
the discrepancy between the binned and estimated broad- 
line AGN BHMFs is as large as two orders of magnitude 
in the mass range around Mbh = lO 9 Af at z = 1.4 
(IShen fc Kellvll2012f) . Therefore, in order to examine the 
cosmological evolution of BHMFs and ERDFs, a sam- 
ple of broad-line AGNs with fainter detection limits is 
needed. 

In order to reveal accretion onto SMBHs in an era of 
violent growth, we examine the BHMF and ERDF of 
broad-line AGNs at z = 1.4 using a sample constructed 
from the X-ray survey of the Subaru XMM-Newton Deep 
Survey (SXDS) field. As sug gested by the SMBH growth 
curves (M arconi "et~alll200i r. a significant part of the ac- 
cretion growth of SMBHs is thought to be taking place 
in the redshift range between 1 and 2. Therefore the 
direct determination of the BHMF and ERDF in this 
redshift range is of critical importance. Thanks to the 
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Redshift 

Figure 1. Redshift distribution of the sample. The thick solid 
line is for the entire sample with spectroscopic or photometric red- 
shifts. AGNs with spectroscopic redshifts are shown with short- 
dashed line. The thin solid line is the broad-line AGNs with spec- 
troscopic redshifts. The dotted line is the broad-line AGNs with a 
black hole mass estimate from either the Mg II or Ha broad-line. 



moderately deep detection limit and wide area of the 
survey, we can construct a large sample of broad-line 
AGNs which is one order of magnitude fainter than that 
available from SDSS. Furthermore, the sample covers the 
flux range around the knee of the X-ray log N-log S re- 
lation (1 x 10~ 14 erg s~ : cnr 2 in the 2-10 keV band; 
ICowie et al.l I2002h . and the sample size is several times 
larger than the deep Chandra surveys in this flux range. 
The sample represents the population of SMBHs that 
dominates the accretion growth of the SMBHs. 

This paper is organized as follows. The details of the 
sample are described in Section [2] Measurements of 
broad-line width of Mg II and Ha lines and Mbh es- 
timates based on these line widths are discussed in Sec- 
tion [3] In Section 21 estimates of the bolometric lumi- 
nosity and Eddington-ratio are presented. In Section O 
we first present the binned BHMF and ERDF calcu- 
lated applying the V max method to the sample of broad- 
line AGNs between z = 1.18 and 1.68 with black hole 
mass estimates. We then present the detection-limit cor- 
rected broad-line AGN BHMF and ERDF derived with 
the Maximum Likelihood method assuming functional 
shapes of the BHMF and ERDF. In this paper, we do not 
include the effect of the uncertainties of the virial black 
hole mass estimate in the BHMF and ERDF determina- 
tion. In Section [SJ the shapes of the corrected broad- line 
AGN BHMF and ERDF are com pared with those in a 
similar redshift range from SDSS dShen et al.ll20ll and 
thos e in the local Universe fro m the ESO/Hamburg sur- 
vey (jSchulze fc Wis otzki 20 ll|). The contribution to the 
binned active BHMF from obscured narrow-line AGNs is 
also discussed. Throughout this paper we adopt the fol- 
lowing cosmological parameters: Ho = 70 km s _1 Mpc -1 , 

fiivr = 0.3, and Qa = 0.7. M agnitu des are given in 

the AB magnitude system (|Oke Ill974l ) unless otherwise 
noted. 

2. SAMPLE 

A sample of z ~ 1.4 AGNs is co nstructed from X - 
ray observations of the SXDS field (jUeda et all 120081 ) . 
The field was observed with the XMM-Newton covering 
one central 30' diameter field at a depth of lOOks expo- 
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Figure 2. Redshift vs. absorption-corrected 2-10 keV luminosity 
of X-ray-selected AGNs in the SXDS. Filled circles and squares rep- 
resent broad-line and narrow-line AGNs with spectroscopic identifi- 
cation, respectively. Broad-line (narrow-line) AGN candidates with 
photometric redshifts are shown with open circles (open squares). 



sure and six flank ing fields with 50ks exposure time each 
(|Ueda et al.H2008D . From summed images of pn, MOS1, 
and MOS2 detectors, there are 866 and 645 sources de- 
tected in the 0.5-2 keV (soft) and 2-10 keV (hard) bands, 
respectively, with a detection likelihood, which is deter- 
mined by point spread function fitting, larger than 7, 
which corresponds to a confidence level of 99.9%. In 
this paper, only X-ray sources in the region covered with 
deep optical imaging data taken with Supr ime-cam on 
the Subaru telescope (|Furusawa et al.ll2008l ) are consid- 
ered. There are 781 and 584 sources in the soft- and 
hard-bands, respectively. Once Galactic stars and clus- 
ters of galaxies candidates are removed, 733 (576) sources 
remain as AGN candidates. Hereafter we call the former 
(latter) sample the soft- (hard-) band sample. The detec- 
tion limit of the survey corresponds to a flux of 6 x 10~ 16 
erg s _1 cm~ 2 (3 x 10~ 15 erg s _1 cm~ 2 ) in the soft (hard) 
band. The area covered with the flux limit is 0.05 deg 2 , 
and 1.0 deg 2 is covered at the flux limit of 4 x 10 -15 erg 
s _1 cm" 2 in the soft band. Considering sources common 
to both samples, there are 896 unique AGN candidates 
in total (see Table [1]). 

In order to identify the X-ray sources spectroscopically, 
optical observations were conducted with various multi- 
object spectrographs on 4m- and 8m-class telescopes. 
The optical spectroscopic observations cover 590 out of 
the 896 sources in total. Even though the observations 
do not cover the entire sample, they are not heavily bi- 
ased toward a specific type of object, since we do not 
use any further discrimination such as color in the tar- 
get selection in most of the observations. Details of the 
observations are summarized in Akiyama et al. (2012, in 
preparation). 

Additional intensive NIR spectroscopic observations 
were made with the Fiber Multi Object Spectrograp h 
(FMOS) on the Subaru telescope (jKimura et al.ll2010l) . 
This instrument can observe up to 200 objects simultane- 
ously over a 30' diameter FoV in the cross-beam switch- 
ing mode with two spectrographs. The spectrographs 
cover the wavelength range between 9000A and 18000A 
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Table 1 

Summary of the SXDS broad-line AGN sample 





N 


N so ft 


Nh ar d 


X-ray sources 


945 


781 


584 


AGN candidates 


896 


733 


576 


Optical spec, observed 


590 


517 


396 


FMOS spec, observed 


851 


704 


548 


Spec, identified 


586 


514 


397 


Mg II broad-linc 


186 


181 


137 


Ha broad-line 


81 


78 


68 


Mg II and Ha broad-line 


52 


51 


11 



Note 

within deep Suprime-cam image coverage 

among the 896 AGN candidates 
among the 896 AGN candidates 

z range 0.489 - 2.329 
z range 0.634 - 1.655 



Note. - 



Table 2 

Summary of SXDS AGNs in the redshift range 1.18 < z < 1.68 





N 


N soft 


Nhard 


Broad-linc AGN sample 


Total (excluding broad-linc AGNs with photometric redshift) ^ 


118 


112 


90 


M BH with Mg II broad-line FWHM 


93 


90 


67 


Additional M B h with Ha broad-line FWHM 


23 


21 


21 


Spcctroscopically-identificd broad-line AGN w/o Mg II or Ha measurement 


2 


1 


2 


Broad-line AGN in the redshift range identified with photometric redshift a 


10 


10 


5 


Narrow-line AGN sample 








Total 


158 


120 


101 


Spcctroscopically-identificd narrow-line AGN 


66 


52 


41 


Narrow-line AGN with photometric redshift 


92 


68 


60 



a Non-detection of broad emission line in the observed wavelength range implies that most of the 
broad-line AGNs with photometric redshift may not be in the redshift range. Therefore we do not 
include them from the total number in the first line and the derivation of BHMF and ERDF, but 
show their number explicitly in the 5th line for reference. See text in Section 2. 



with spectral resolution of R ~ 800 at A ~ 1.55/zm in 
the low-resolution mode. A total of 851 sources were 
observed with this setup during guaranteed, engineering 
and open-use (S11B-098 Silverman et al. and S11B-048 
Yabe et al.) observations. The optical and NIR ob- 
servations spectroscopically identify 586 out of the 896 
sources. The optical and NIR spectra obtained in the 
identification observations are used for the broad-line 
width measurements described in the next section. 

The remaining 310 sources cannot be identified spec- 
troscopically, mostly because of their faintness. Most of 
them are fainter than i = 23.5 magnitude. For such ob- 
jects, photometric redshifts have been estimated using 
the H yperZ photometric redshift code dBolzonella et al.l 
120001 ) with galaxy and QSO Spectral Energy Distribu- 
tion (SED) templates. Photometric data in 15 bands 
covering from 1500A to 8.0/xm are used in the estima- 
tion. In order to reduce the number of AGNs with signif- 
icantly different photometric redshift from spectroscopic 
one ("outliers"), we apply two additional constraints in 
addition to the % 2 minimization considering the prop- 
erties of the spectroscopically-idcntificd AGNs. First 
one is that the objects with stellar morphology in the 
deep optical images are z > 1 broad-linc AGNs. Al- 
most all X-ray sources with stellar morphology are iden- 
tified with broad-line AGNs at z > 1 in SXDS. They 
show a bright nucleus and their observed optical light 
is dominated by the nuclear component. Second one is 
the absolute magnitude range of the galaxy and QSO 



templates. Considering the absolute magnitude range 
of spcctroscopically-identificd broad-line and narrow-line 
AGNs, we limit the z-band absolute magnitude range of 
the galaxy (QSO) template between M z = -20.0- -25.0 
(mag) [M z = -22.0- -26.5 (mag)). 

The accuracy of the photometric redshifts are exam- 
ined by comparing them with the spectroscopic redshifts. 
The median of Az/(1 + z spcc ) is 0.06 for the entire sam- 
ple. We further examine the accuracy by the normal- 
ized median ab s olute deviation (NMAD; a z ) following 
IBrammer et alj (|2008f ). For the entire sample, a z is 
0.104, which is larger than that of the photometric red- 
shift estimati ons for X-ray-selected AGNs with medium 
band filters (|Cardamone et all [2010l ; iLuo et al.l [20100 . 
The a z for broad-line AGNs (0.201) is larger than that 
for narrow-line AGNs (0.095). This is because there is 
no strong feature in the SEDs of the broad-line AGNs 
except for the break below Lya. 

From the photometric redshift determination, not only 
their photometric redshifts, but also their types of SED 
can be constrained. For spectroscopically identified 
AGNs, there is a good correlation between the spectral 
type and the best-fit SED type; narrow-line and broad- 
line AGNs are fitted well with galaxy and QSO tem- 
plates, respectively Therefore, we classify objects fit- 
ted better with the QSO templates as broad-line AGNs 
and the others as narrow-line AGNs. The classifica- 
tion does not perfectly match the spectroscopic classi- 
fication: in the redshift range, 10 out of 66 (31 out 
of 118) spectroscopically-identified narrow- line (broad- 



Broad-line AGN BHMF and ERDF at z=1.4 



5 



line) AGNs are photometrically classified as broad-line 
(narrow-line) AGNs. The SEDs of the spectroscopically- 
unidentificd objects suggest that most of them are ob- 
scured narrow- line AGNs above redshift 1. For 6 objects, 
no photometric redshift can be estimated because they 
are detected only in a few bands. They are faint and 
are unlikely to be broad-line AGNs in the redshift range 
between 1.18 and 1.68. Further details of the photomet- 
ric redshift determination is discussed in Akiyama et al. 
(2012, in preparation). 

In this paper, for objects with spectroscopic identifi- 
cation, we designate objects as broad-line AGNs if they 
show either Mg IIAA2796, 2803 or Ha emission lines with 
width greater than 1000 km s _1 . We estimate black hole 
mass of the broad-line AGNs with either broad Mg II or 
Ha line. The threshold is narrower than typical thresh- 
old used to discriminate broad-line AGNs (1500 or 2000 
km s" 1 ). We determine the threshold considering the 
distribution of the FWHM of the broad-line AGNs i n 
the local universe (jHao et al.ll2005t iStern fc Laorll2012ft . 
The broad-line AGNs with the line FWHM close to 
the threshold correspond to the narrow-line Seyfert Is. 
Broad Mg II and Ha lines are detected for 186 and 81 
AGNs respectively, with redshifts in the range between 
0.5 and 2.3. For 52 objects, both broad Mg II and Ha 
lines are detected. For 29 objects, a broad-line is only 
detected in Ha. This is mostly due to the lack of op- 
tical spectra. Only 4 AGNs (SXDS0215, SXDS0387, 
SXDS0527, SXDS0728) with broad Ha line show no 
broad Mg II line, although their optical spectra cover 
the Mg II wavelength region and are deep enough to de- 
tect continuum emission. Such AGNs are thought to be 
moderately affected by dust extinction and we correct for 
this in the determination of their continuum luminosity 
for the black hole mass and the intrinsic bolometric lu- 
minosity estimates. Details are given in Section [3] and 
® 

For the derivation of the broad-line AGN BHMF, we 
limit the sample to the redshift range between 1.18 and 
1.68. The FMOS i7-band observation covers the broad 
Ha line for AGNs in the redshift range. Considering 
the typical detection limit for the broad Ha line in 
this redshift range (£h q = 2 x 10 42 (erg s -1 )), we ex- 
pect a broad Ha can be detected for broad-line AGNs 
brighter than L 2 -io kcV = 3 x 10 43 (erg s _1 ), if they 
have the Ha to h ard X-ray lumino sity ratio typical of 
broad-line AGNs (|Ward et al.1 11988ft . Photometric red- 
shift estimates suggest 10 of the unidentified sources can 
be broad-line AGNs in this redshift range. All but one 
of the 10 candidates have estimated hard X-ray lumi- 
nosity larger than 3 x 10 43 (erg s _1 ). But, 8 of the 
10 broad-line AGN candidates do not show a broad- 
line in the FMOS observations. Considering the un- 
certainty of photometric redshift for broad-line AGNs, 
we expect they are broad-line AGNs at outside of the 
redshift range. However, th ere is a 0.4 dex sc atter be- 
tween the L2-10 kcV and Lw^ tjWard et al. ill 988ft . and it is 
still possible that they have weaker broad Ha line than 
the typical broad-line AGNs. Furthermore, additional 
one broad-line AGN candidate does not have broad-line 
in the optical spectroscopy, it may also be a broad-line 
AGN at outside of the redshift range. Considering the 
non-detection of broad-line in the observed wavelength 



range, we do not include the 10 photometric candidates 
of broad-line AGNs in the redshift range in the deriva- 
tion of the BHMF and ERDF below. The numbers of 
X-ray-selected AGNs in the redshift range are summa- 
rized in Table [2l The median redshift of the sample is 
1.43. 

The redshift distribution of the sample is shown in Fig- 
ure [T] for the redshift range 0.5 to 2.5. The thick solid 
line shows the redshift distribution of the X-ray AGNs 
with spectroscopic or photometric redshifts. The dashed 
line shows the distribution of spectroscopically-identified 
AGNs. The thin solid line shows the distribution for 
broad-line AGNs with spectroscopic redshifts. The dot- 
ted line is the distribution for broad-linc AGNs that have 
black hole mass estimates with either broad Mg II or Ha 
emission lines. 

In Figure [2j the absorption-corrected 2-10 keV lumi- 
nosity of AGNs are show n as a function of redshift. Fol- 
lowing |Uedi~eFII] (|2003ft , we assume the intrinsic shape 
of the X-ray spectrum and estimate the intrinsic column 
density from the observed hardness ratio and redshift. 
The intrinsic X-ray spectrum of AGNs is modeled by a 
combination of a power-law component with high-energy 
cut off E~ T x exp (—E/E c ) and a reflection component. 
A photon index T of 1.9 and cutoff energy E c of 300 keV 
are assumed. We calculate the reflection com ponent with 
the "pexrav" (Magd ziarz fc Zdziarski l995!) model in the 
XSPEC package assuming a solid angle of 27r, an incli- 
nation angle of cosi = 0.5, and solar abundance of all 
elements. Strength of the reflection component is about 
10% of the direct component just below 7.1 keV. The 
intrinsic SEDs are modified with intrinsic photo-electric 
absorption described by a hydrogen column density, TVh . 
The JVh value of each object is evaluated from the ob- 
served 0.5-2 keV and 2-4.5 keV hardness ratio. The 
absorption-corrected luminosity in the 2-10 keV band is 
derived from the observed 2-10 keV count-rate by cor- 
recting for the photo-electric absorption. For objects 
only detected in the 0.5-2 keV band, the count-rate in 
that band is used instead of the 2-10 keV count-rate. 
Most of the broad-line AGNs have a hardness ratio con- 
sistent with no significant absorption and the required 
correction is small. Only 5 out of the 215 broad-line 
AGNs have logA^ as large as 23; for these the correc- 
tion in luminosity is as large as 0.5 dex. 

The absorption-corrected 2-10 keV luminosity of 
AGNs distribute between L 2 -io koV = 10 43 and 10 45 (erg 
s- 1 ) in the redshift range of z = 1 - 2. The SXDS AGNs 
cover the most important part of the accretion growth 
of the SMBHs. AGNs in the luminosity range domi- 
nate the hard X-ray luminosity density of the Universe 
in the redshift range, furthermore, the hard X-ray lumi- 
nosity density as a function of redshift peaks at z = 1 
(|Aird et al.ll2010ft . 

3. LINE WIDTH AND LUMINOSITY MEASUREMENT 

3.1. Method for Black Hole Mass Estimation 

Assuming that the scaling relationship between lumi- 
nosity and broad-line region size derived by reverbera- 
tion mapping for local broad-line AGNs is applicable to 
broad-line AGNs at high-redshifts, black hole masses of 
broad-linc AGNs can be estimated from their continuum 
luminosities and line widths of the broad lines. We es- 
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timate the black hole mass of z ~ 1.4 broad-line AGNs 
with Mg II broad-line widths measured in the optical 
spectra. As the optical spectroscopy does not cover the 
entire sample, we supplement these with Ha broad-linc 
widths measured in the NIR spectra. In this subsection, 
we first introduce the equation used to estimate the black 
hole mass. 

There are several calibrations available for the 
black hole mass estimat io n usin g the Mg II broad 
line (IMcLure fe Jarvisl l2002t IMcLure fe Dunk) ' 
2001 IMcGill et all I2008L IVestergaard k Osmerl 1200 



Wang et alJI2009t IShen et al.H2011HRafiee k Halll[201ll) 



We use the black hole mass estimate from the Mg II 
broad-line FWH M calibrated by broad - line QSOs in 
the SDSS DR3 (IVestergaard fe Osmerl [2009h . The 
estimation is consistent to within 0.1 dex of the H/3 and 
C IV mass estimates with single-epoch spectra. They are 
calibrated to the black hole mass from th e reverberation 
mapping (jVestergaard k Peterson! [2006) . It needs to 
be noted that the black hole mass determined with 
the single-epoch H/3 spectrum typically has scatter of 
0.4-0.5 dex around that from t he reverberati o n map ping 
(IVestergaard fe Peterson[|2006f) . IPark et all (|2012f ) also 
estimate the uncertainty of the single-epoch black hole 
mass to be 0.4-0.5 dex. Furthermore, if AGNs are 
close to the Eddington limit, the influence of radiation 
pressure may cau se an underestimatio n of the black hole 
mass by 0.5 dex ([Marconi et al.ll2008D . 

Because the 3000A monochromatic luminosity is avail- 
able for most of the broad-linc AGNs with spectro- 
scopic observations, we use the following equation from 
IVestergaard k Osmerl ([20091 ) incorporating the 3000A 
monochromatic luminosity, £a30oo- 



m bh [m ; 
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6.86 



FWHM 



Mgll 



1000 km s" 



A3000L 



A3000 



10 44 erg s" 



In this calibration, the FWHM of broad Mg II line, 
FWHMMgii, is measured by fitti ng multiple Gaussians 
to the Mg II emission line profile. IVestergaard k Osmerl 
(|2009t ) remove narrow-line com ponent of Mg II if nec- 
essary ( Vestergaard et ahl 120111) . In other calibrations, 
such as lMcLure fe Jarvisl ([20021 ) and IMcLure fc Dunlopl 



(|2004h . fit the Mg II profil e with a single broad -line and 
narrow-line components. iRafiee fc~Tl all (|2011|) use the 
line dispersion, a, of the broad- line because a correlates 
more tightly with the delay observed in reverberation 
mapping , i.e. size of the broad -line region, than with the 
FWHM ([Peterson et al.ll200l . 

In Equation (1), the broad-linc region size R is as- 
sumed to follow L a with a of 0.5, which is equivalent 
to assuming that the broad-line regions of various AGNs 
can be described as having similar ionization states, ion- 
izing photon spectra and electron densities. The value 
of a, based on reverberation mapping results for H/3 
broad-lines, i s determined to be 0-47, 0.62 ± 0.14, and 
0.518 ±0.039 (IMcLure fc: Jarvisl[2002l: I McLure fc Dunlod 
120041 : iBentz et al.ll2006l) . respectively, all consistent with 
a of 0.5 within the uncertainties. 

Regarding the coefficient for the virial product (e of 
M BH = eFWHM H/ 3 2 £ ' 5 or / of M BH = f^npL - 5 ), 
Equation (1) is based on the calibration done by 



lOnken et"aT1 ((20041) (e of 1.4 or / of 5.5) assuming the 
estimated Mbh of local broad-line AGNs from reverber- 
ation mapping and their bulge velocity dispersions fol- 
low the black hole mass and the bulge velocity disper- 
sion relation of non-active galaxies in the local Universe. 
The coefficient is consistent wit h that derived w ith lo- 
cal Seyfert 1 galaxies (/ = 5.2; IWoo et alll20To[) . Re- 
cent calibration shows broad-line AGNs hosted in barred 
galaxies a re consistent with si gnificantly smaller values 
(/ = 2.3: IGraham"et~aT1 l20ll . but the range of black 
hole mass of the barred galaxies is smaller than the cur- 
rent sample, and non-bar red galaxies with larg er Mbh is 
consistent with / ~ 5.4 (|Graham et al.l [20111 ). thus we 
use Equation (1). 

The optical spectroscopic observations do not cover all 
of the broad-line AGNs at z = 1.18 - 1.68; 25 out of 
118 spectroscopically- identified broad-line AGNs in the 
rcdshift range do not have Mg II data. Therefore, we 
also utilize the Ha FWHM in addition to the Mg II 
FWHM; additional 23 broad-line AGNs have Ha broad- 
line data. Although the FMOS spectra cover H/3 in J- 
band, the strength of the H/3 broad-line is 3 times or more 
weaker than the Ha broad-line and the uncertainty of the 
FWHM of the broad H/3 is significantly larger than that 
of the broad Ha line. Therefore, we do not use the H/3 
FWHM. Because ionization potentials of hydrogen and 
Mg II are similar, Balmer and Mg IIAA2796, 2803 broad 
lines are expected to be emitted in a similar region. A 
detailed photoionization model calculation indicates that 
the equivalent widths of Ha and Mg II lines have a sim- 
ilar dependency on the cloud density and ionization pa- 
ramete r, i.e., they are em itted from similar broad-line 
clouds (jKorista et al.lll997h . Considering this similarity, 
we use the same black hole mass equation employed for 
Mg II FWHM above for Ha FWHM, after correcting for 
a small systematic difference between Mg II FWHM and 
Ha FWHM, as detailed below. We do no t use the scal- 
ing r elation calibrated for Ha broad-line (|Greene fc Hoi 
120051) in order to be consistent within our sample. The 
derivation of the 3000A monochromatic luminosity is dis- 
cussed in Section [331 



3.2. Mg II Line Width Measurements with Optical 
Spectrum 

Optical spectra of the AGNs were obtained with vari- 
ous instruments on 4-8m class telescopes such as 2dF on 
the Anglo-Australian Telescope, VIMOS and FORS on 
the Very Large Telescope, FOCAS on the Subaru tele- 
scope, DEIMOS on the Keck telescope, and IMACS on 
the Magellan telescope. Most of them were obtained 
with spectral resolution of R ~ 250 — 500. Details of 
the observations are described in Akiyama et al. (2012, 
in preparation). The optical spectroscopic data were re- 
duced using standard procedures and are corrected for 
the dependence of the sensitivity on wavelength with 
standard star observations. We further correct the nor- 
malization of the spectra to match the observed i?-band 
magnitudes in the deep Suprime-cam images. The opti- 
cal photometric data are corrected for the Galactic ex- 
tinction in the SXDS field (A R of 0.07 mag). We do not 
correct the optical spectroscopic data for the wavelength 
dependence of the Galactic extinction which is negligibly 
small. This correction docs not affect the line width mca- 
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Figure 3. Examples of broad Mg II fitting. Upper panels show the observed data (thin solid line) and the best fit model from the 
power-law continuum and Fe II fitting (thick solid line). Power-law component of the best fit model is also shown separately. Middle panels 
show the results of the broad Mg II line fitting. Pure Mg II component after subtracting the power-law continuum and Fe II components 
and best fit model from the Mg II fitting are shown with thin solid line and thick solid line, respectively. Each component of the best fit 
model is also shown. Only the wavelength range used for the broad-line fitting is plotted. Bottom panels show the residual after fitting. 
Left) SXDS0814 with single broad-line, and right) SXDS0749 with 2 components. 
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Figure 4. FWHMs of the broad Mg II emission line determined 
after Fe II fitting with with VW01 and T06 Fe II templates. The 
dashed line shows the equal FWHM line. 

surement, but affects the continuum flux measurement. 
The normalization can be affected by the variability of 
broad-line AGNs between the epochs of the imaging and 
the spectroscopic observations. Typically there is one 
year gap between the imaging and spectroscopic obser- 
vations . The structure func tion of optical variability of 
QSOs (jCristiani et al.lfl996f) suggests that there can be 
0.2 mag variation during the time lag on average. There- 
fore, the uncertainty of Mbh due to the time-variation 
of AGNs is expected to be 0.04 dex. 

In order to determine the Mg II FWHMs of broad-line 
AGNs, it is necessary to consider Fe II emission lines 
as well as the power-law continuum component in the 
UV wavelength range. Because there are many broad 
Fe II emission lines in the wavelength range, they look 
like an additional continuum component to the power- 
law continuum of broad-line AGNs. We use a Fe II tem- 
plate derived from the U V spectrum of the narrow-lin e 
Seyfert 1 galaxy, I Zw 1 (jVestergaard fe Wilkes! 120011 ). 



This template covers the rest-frame wavelength range be- 
tween 1074A and 3089A. We do not include the Balmer 
continuum in the fitting, because the wavelength cover- 
age is not wide enough to constrain its contribution. The 
ignorance of the Balmer continuum does not affect the 
Mg II FWHM measurements significantly, but the lumi- 
nosity of the power-law continu um can be overestimated 
by 0.12 dex (|Shen fe Lid 120121) . 

A fit to the power-law continuum, Fe II emission lines, 
and Mg II emission line is carried out as follows. First, 
we determine the normalization of Fe II and continuum 
component using x 2 minimization in the two rest-frame 
wavelength ranges, 2500-2700A and 2900-3000A in which 
Mg II broad-line component is negligible. These are 
close to the pure Fe emissio n windows nos.9 and 10 in 
iVestergaard fe Wilkes! (|2001f ). We vary the line width 
of the Fe II emission lines from 1000 km s _1 to 15000 
km s _1 with a step size of 250 km s _1 by convolving 
a Gaussian profile with the Fe II template which has a 
velocity width of 900 km s _1 . We assume a constant 
line width for all Fe II emission lines. The scaling of the 
Fe II emission is changed from to 100% of the observed 
continuum level with a step size of 0.01%. The contin- 
uum component is modeled with a power-law spectrum 
(/„ oc v~ a ). The observed wavelength ranges affected by 
strong night sky lines (5555-5605A and 6270-6320A) are 
removed in the fitting. Because some optical spectra do 
not have noise level estimation from the standard reduc- 
tion method, we estimate the noise level for each spec- 
trum as follows. Considering that the noise level does not 
vary significantly within the observed wavelength range, 
we use a constant noise level for the entire wavelength 
range. In the first stage of the fitting, we use the stan- 
dard deviation determined within the wavelength range 
as the noise level. Subsequently, we determine the noise 
level from the rms of the residual of the first fitting, and 
carry out the final fit for the continuum. The noise level 
is used for the Mg II line profile fitting as well. Exam- 
ples of continuum fits are shown in the upper panels of 
Figure [3l 

By subtracting the Fe II emission and power-law con- 
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tinuum components, the Mg II broad-linc component 
is extracted. Then we measure the FWHM of the 
Mg II broad-emission line after fitting its line profile 
with multiple Gaussians using \ 2 minimization. In the 
fitting, we do not consider the doublet component of 
Mg IIAA2796, 2803, because the separation is small and 
does not affect the measured width of the broad Mg II 
line. We use mpfit packag e for python for \ 2 min- 
imization (|Markwardtl [20091. Mpfit uses Levenberg- 
marquardt algorithm to derive the best fit parameters. 
For some objects, the Mg II line profile cannot be de- 
scribed with a single Gaussian component. In such cases, 
we consider up to three broad Gaussians for the broad 
line. If necessary, we include narrow doublet absorption 
lines. Sometimes we also include a narrow Gaussian com- 
ponent to remove artificial spiky noise features. Because 
no object shows a significant existence of narrow Mg II 
line, we do not include a narrow-emission line component 
in the fitting. No inclusion of the narrow-line component 
differs from th e fitting method u s ed in most of the liter- 
atures such as iMcLure fc Jarvisl (|2002[ ) . Once the pure 
Mg II broad-line component is fitted with the multiple 
Gaussian components, the FWHM of the broad Mg II 
line is measured with the best fit profile constructed by 
combining the multiple Gaussians. We do not include 
absorption lines in the combination. We introduce mul- 
tiple Gaussian components in order to reproduce the ob- 
served Mg II broad-line profile smoothly and here we are 
not concerned with the physical meaning of the difference 
from the single Gaussian profile. Examples of the result- 
ing Mg II broad-line fitting are shown in the middle and 
bottom panels of Figure [3J Because the spectra are ob- 
tained with various instruments, the spectral resolution 
of data varies from object to object. The spectral resolu- 
tion is evaluated for each spectrum by using line width of 
the arc lamp spectrum or night sky emission lines. The 
measured FWHMs are corrected for the intrinsic spectral 
resolution. 

The uncertainty of the FWHM values for the combined 
multiple Gaussian profile is not available from the fit- 
ting with mpfit package (uncertainty is only available 
for each Gaussian component). Therefore we evaluate 
the uncertainty for the FWHM of each object from the 
rms scatter of FWHMs measured in mock spectra con- 
structed from the best fit profile of the object. We con- 
struct the mock spectra as follows. First, the best fit 
multi-Gaussian model is shifted by several pixels in wave- 
length from its original position. Then the shifted model 
profile is combined with the residual of the original fit- 
ting. By monotonically increasing the shift and changing 
the sign of residual in each pixel randomly, we construct 
100 mock spectra. The mock spectra are fitted in the 
same way as the original data and FWHMs are mea- 
sured. Finally the rms scatter of the derived FWHMs is 
used as the uncertainties of the FWHM measurement. 
For AGNs whose Mg II broad-line is fitted with sin- 
gle Gaussian component, we compare the uncertainties 
derived from the \ 2 statistics and from the scatter of 
the mock measurements. They are consistent with each 
other, although the rms scatter of the mock measure- 
ments is slightly smaller than the uncertainties from the 
X 2 statistics. Hereafter, we use the uncertainty derived 
with the scatter of the mock measurements. The fitting 
results are summarized in Table |3l Column 4 of the ta- 



ble describes the model used for each object; "OncBL", 
"TwoBL" and "ThrccBL" indicate fitted with one, two 
and three broad lines respectively. "OncBLOncAbs" in- 
dicates a model with one broad and one absorption line. 

For the Mg II profile fitting, we also use the specfit soft- 
ware in the stsdas package of IRAF. This package uses 
Marquardt algorithm or simplex algorithm for \ 2 mini- 
mization. We compare the results obtained with mpfit 
and specfit for each object. The rms scatter of the dif- 
ference in FWHMs from the two measurements is 0.06 
dex. The resulting uncertainty in -Mbh due to the scat- 
ter is 0.12 dex. We use the results obtained with mpfit 
hereafter. 

The measured FWHM of Mg II broad-line can be af- 
fected by the template used for the Fe II fitting. For the 
Fe II fitting, we use the empir i cally d erived Fe II template 
from iVestergaard fc Wilkes! (pOOll ) (VW01 template). 
There is no Fe II emission in the wavelength range be- 
tween 2770A and 2820A in the template. iTsuzuki et"aT1 
(|20061 ) also derived the Fe II template (T06 template) 
from UV and optical spectra of I Zw 1. Utilizing the 
wide wavelength coverage available, they fit the contin- 
uum with a power-law and Balmer continuum and also 
utilize the Ha line profile for removing the Mg II line 
component. The important difference between the two 
templates is the Fe II contribution to the blue wing of 
the broad Mg II line. In the T06 template, the excess 
wing seen on the blue side of the Mg II line of I Zw 1 
compared with its Ha line profile is regarded as a contri- 
bution from the Fe II emission line at around 2790A. In 
order to examine the effect of different Fe II templates on 
the FWHM measurement, we also apply the Mg II emis- 
sion line fitting process described above using the T06 
template for 95 objects whose broad Mg II lines are fitted 
well with one broad-line component with the VW01 tem- 
plate. Figure S] compares the log FWHMs values derived 
with the two templates. There is a systematic offset of 
0.04 dex; the FWHMs derived with the T06 template are 
systematically smaller than those with the VW01 tem- 
plate. The rms scatter of the difference is 0.05 dex after 
removing the systematic offset. The resulting systematic 
uncertainty of Mbh is 0.08 dex. In order to compare the 
broad-line AGN BHMFs from literature, we follow the 
same fitting procedure using the VW01 template, but it 
needs to be noted that the FWHMs can have a system- 
atic uncertainty due to the Fe II template difference. 

3.3. Ha Line Width Measurements from NIR Spectrum 

The NIR spectra of X-ray sources obtained with FMOS 
were reduced with the pipeline da ta reduction software, 
FIBRE-pac (jlwamuro et all l2012f ). The resulting 1-d 
spectra are corrected for atmospheric absorption and sen- 
sitivity dependence on wavelength using relatively bright 
F-G type stars observed simultaneously with the targets. 
Because we do not use the absolute flux of the continuum 
component, we do not apply further correction for nor- 
malization based on the photometry. Noise level of each 
spectrum is also estimated in the pipeline data reduction 
software. 

The profile of the broad Ha line is fitted with mpfit 
in a similar way for broad Mg II component. The pos- 
sible existence of the narrow emission lines of Ha and 
[N II]AA6583,6548 makes the fitting more complicated 
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Table 3 

Broad-line fitting results for the broad-line AGNs 



Mg II 



Ha 



ID 


z 


log[-^HX 


FWHM 


model 


FWHM 


model 


log[A3000L A30 oo 


R-i 






/ [erg s" 1 ] ] 


[km s" 1 ] 




[km s _1 ] 




/ [erg s" 1 ] ] 


[mag] 


0010 


1.225 


44.05 


5341 ± 12 


TwoBL 


5135 ± 116 


OneBL 


44.89 


-0.04 


0018 


1.452 


44.42 






5488 ± 92 


BLNL 


45.16 


-0.02 


0019 


1.447 


44.66 


4823 ± 312 


OneBL 






45.41 


-0.06 


0023 


1.534 


44.14 






5602 ±112 


OneBL 


45.10 


0.48 


0027 


2.067 


43.75 


4136 ± 1499 


TwoBL 






45.85 


0.27 


0034 


0.952 


43.55 


4286 ± 1005 


TwoBL 


2459 ±123 


TwoBL 


44.72 


-0.15 


0036 


0.884 


44.16 


3326 ± 22 


TwoBL 


2790 ± 50 


TwoBL 


45.24 


0.25 


0037 


1.202 


43.81 


4516 ± 70 


TwoBL 






44.45 


-0.06 


0050 


1.411 


44.03 


2046 ± 120 


OneBL 


1800 ± 68 


OneBL 


44.80 


0.14 


0056 


1.260 


44.19 






8171 ± 456 


BLNL 


44.74 


-0.03 



a Table 3 is published in its entirety in the electronic edition of ApJ. A portion is shown here for guidance 
regarding its form and content. 
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Figure 5. Examples of broad Ha fitting. Upper panels show the observed data (thin solid line) and the best fit model with each 
component (thick solid lines). Lower panels show the residual from the fitting (thin solid line). Thick solid lines in the panel enclose the 
estimated la noise level at each wavelength. Left) SXDS1097 with single broad-line, middle) SXDS0817 with 2 broad-lines, and right) 
SXDS0018 with broad-line and narrow Ha and [N II] AA6548, 6583 lines. 



than for Mg II broad line. We include the narrow emis- 
sion lines in the fitting only if a prominent narrow emis- 
sion feature, such as asymmetric profile due to the 1:3 
flux ratio of [N II] lines, is observed. For the broad Ha 
component, we use up to two Gaussians to fit their pro- 
file. In the fitting process, we assume a constant contin- 
uum component because the continuum of the observed 
NIR spectra does not show significant tilt in the region 
around Ha emission line. Due to the existence of the OH 
suppression mask, the strength of the narrow emission 
line can be underestimated even after the sensitivity cor- 
rection process if the narrow emission line is close to the 
masked wavelength. We apply the fitting procedure tak- 
ing into account the underestimation due to the optical 
masking and the sensitivity correction proce ss. The de- 
tails o f the fitting procedure are described in lYabe et al.l 
(|2012f ). It needs to be noted that the effect of OH sup- 
pression mask only affects the narrow-emission line and 
the effects on the broad-emission lines and continuum are 
negligible. Examples of Ha fitting are shown in Figure [5j 
The uncertainty in the FWHM measurements of the 
broad Ha emission line is evaluated in the same way as 
for the broad Mg II emission lines. We construct 10 



model spectra by adding shifted best fit multi-Gaussian 
model to the residual of the fitting. Then, the same fit- 
ting procedure is applied to the model spectra and the 
rms scatter of the measured FWHMs is used as the un- 
certainty of the FWHM measurement. In this fitting 
process we also take into account the effect of the OH 
suppression mask. 

The resulting FWHM of the broad Ha emission is com- 
pared with that of the broad Mg II line for broad-line 
AGNs with FWHM measurements for both lines in Fig- 
ure El We remove Mg II FWHMs of 4 broad-line AGNs 
(SXDS0590, SXDS0630, SXDS0790, and SXDS0969), 
because the signal-to-noise ratios of their Mg II spec- 
tra are much lower than those for their Ha spectra. 
Therefore, 48 broad-line AGNs are plotted in the fig- 
ure. The sample is divided with the absorption-corrected 
X-ray luminosity. Broad-line AGNs that are brighter 
(fainter) than L2-iokcV = 10 44 erg s _1 are marked with 
filled squares (crosses). The measured FWHMs of the 
broad Mg II and Ha lines roughly follow the equality 
line. However, there may be a tendency that broad-line 
AGNs with lower luminosity have systematically smaller 
broad Ha FWHM than broad Mg II line. The distribu- 
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Figure 6. Comparison between FWHMs of Ha and Mg II broad- 
lines for SXDS AGNs. Filled circles and crosses represent broad- 
line AGNs brighter and fainter than -L2-10 keV = 10 44 erg s — 1 , 
respectively. The thick solid line represents the BCES bisector 
fitting result for the SXDS AGNs. Open squares, triangles, and 
circles are AGNs from Shen & Liu (2012), Greene et al. (2010b), 
and McGill et al. (2008), respectively. The relationships between 
the Mg II and FWHMs from McLurc et al. (2002) (short 
dashed line), Onken & Kollmeier (2008) (thin solid line), Wang et 
al. (2009) (dotted line), and Croom et al. (2011) (long dashed line) 
are shown after converting H/3 FWHM to Ha FWHM. 



tion is broadly consistent with those of broad-line AGNs 
measured in the literature. We show the M r II and Ha 
FWHMs of individua l broad- line A GNs fromlShen fc Liu! 
(I20T1 . [Greene et all (12010b!) and IMcGill et all (120081) 
The sa mples of iShen fc Liul (|2012D and iGreene et al.l 
(l20T0bl) are luminous broad-line AGNs with Lboi of 
10 46 ~ 48 erg s -1 at z = 1 — 2. On the contrary, the 
IMcGill et all (|2008f ) sample covers broad-linc AGNs with 
Lboi of around 10 45 erg s _1 at z ~ 0.3. The distribution 
of the former samples are consistent each other and follow 
similar trend of the SXDS luminous broad-line AGNs. 
The latter sample has systematic offset from the former 
samples and shows smaller Ha broad-line FWHM than 
Mg II FWHM. The trend is similar to that seen in the 
SXDS less-luminous broad-line AGNs. 

We determine the relationship between the Mg II and 
Ha FWHMs applying a BCES (Bivariate Correlated Er- 
ror and intrinsic Scatter) bisector regression analysis 
(jAkritas fc Bershadvl [1999) to the 48 broad- line AGNs 
including both luminous and less-luminous AGNs. Con- 
sidering the size of the sample, we do not divide the 
sample by luminosity in the analysis. The resulting rela- 
tionship is 

log FWHM Mg ii 

= (0.795 ± 0.075) x logFWHM HQ + (0.771 ± 0.273). 

(2) 

The rms scatter of log FWHMMgii determined with the 
above relationship using the measured log FWHMMgii is 
0.11 dex with a resulting uncertainty for logMeH of 0.22 
dex. The less luminous broad-line AGNs show systematic 
offset of 0.1 dex on average from the relation. 




log[ Corrected L 2 . 10 kev / [erg s" 1 ] ] 

Figure 7. Absorption corrected 2-10 keV luminosity vs. 3000A 
monochromatic luminosity of the SXDS broad-line AGNs (filled 
circles). Broad-linc AGNs redder than R — i > 0.3 magnitude are 
marked with large open circles. The thick solid line indicates the 
relationship determined by Marconi et al. (2004). The dotted line 
shows the relationship determined with the BCES bisector analysis 
for broad-line AGNs bluer than R — i < 0.3 magnitude. 



The relationship betw een the FWHMs of the H/3 and 
Mg II broad -lines from iMcLure fc Jarvisl ()2002[ ) {short 
dashed line) , 
IWang et ail 



Onken fc Kollmeierl ([2008D ( thi n solid line), 
2009P ) (dotted line). and lCrooml (|2011l ) (long 
dashed line) are also shown in the figure. These re- 
lationships for the H/3 FWHMs are converted to those 
for the Ha FWHMs using the relationship between the 
FWHMs of the broad Ha and H/3 lines determined from 
z < 0.35 for local bro ad-line AGNs of SDSS (Eq.(3) 
of IGreene fc Hoi [2005T ) . The relationship between the 
FWHMs of Ha and H/3 is consiste nt with a recent de- 
termination of IShen et ahl ( 2008bfl but is sligh t ly off - 
set from the relationship of ISchulze fc Wisotzkil (|2010D . 
The distribution of SXDS broad-line AGNs inconsis- 
tent with the relation ship o f lMcLure fc Jarvisl (l2002|) and 
lOnken fc Kollmeierl (|2008D . The relationship of |Croom| 
( 20111) has a steeper slope than the other relationships 
and follows the distribution of SXDS broad-line AGNs 
except for th e objects with the l argest FWHMs. The re- 
lationship of I Wang et al.l (j2009f ) shows a systematic off- 
set from the other relationships and the distribution of 
SXDS broad-line AGNs. T he origin of t he sh ift is un- 
clear, but it is possible that IWang et all (|2009l) use T06 
Fe II template, and FWHM of Mg II is measured sys- 
tematically smaller than other measurements with VW01 
Fe II template (see Section l3~2|) . 

We use the relationship shown as the Equation (2) to 
convert the Ha FWHM to Mg II FWHM and the same 
black hole mass equation for Mg II FWHM is applied. In 
total, we use the Ha FWHM measurements for 23 out of 
116 broad-line AGNs at redshifts between 1.18 and 1.68. 

3.4. 3000A Monochromatic Luminosities 

For most of the objects, we derive 3000A monochro- 
matic luminosity using the best-fit power-law continuum 
component described in Section |3T2| We do not include 
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the contribution of the Balmer continuum in the fitting, 
the 3000A monochr omatic luminosity can be overesti- 
mated by 0.12 dex (jShen fc Liull2012h . Optical spectra 
covering rest-frame 3000A are not available for broad- 
line AGNs that only covered by the NIR spectroscopic 
data. We estimate their 3000A monochromatic luminos- 
ity by interpolating multi-band photometry data. All of 
the broad-line AGNs are detected in the deep multi-band 
images obtained with Suprime-Cam. We derive their 
rest-frame 3000A flux by interpolating the photometric 
measurements in the neighboring two bands around rest- 
frame 3000A. The photometric data can include broad- 
emission lines and Balmer continuum as well, and the 
3000A luminosity may thus be affected by the broad-line 
component. In order to estimate this effect, we compare 
the 3000A luminosity derived from normalized spectra 
and multi-band photometry for objects with both mea- 
surements. They are consistent each other within the 
rms of 0.10 dex and resulting uncertainty of M B h is 0.05 
dex. It needs to be noted that the optical spectra are 
normalized to match the i?-band photometry from the 
imaging observations as described in Section [3~2l there- 
fore the scatter only reflects the object-to-object varia- 
tion of the strength of the broad-line components. For 
the estimation of the 3000A monochromatic luminosity 
with the photometric data, the contribution from the 
Balmer continuum is not considered. Uncertainty asso- 
ciated with the variability of broad-line AGNs is already 
described in Section l3~2l 

For mildly obscured broad-line AGNs, the 3000A lu- 
minosity can be affected by dust extinction. Addition- 
ally for low-luminosity broad-line AGNs, the 3000A lu- 
minosity can be affected by a host galaxy component. 
In such cases, the neither the 3000A luminosity derived 
from the spectra nor the photometry is a good indica- 
tor of the intrinsic UV luminosity. In Figure [JJ the 
monochromatic 3000A luminosity and absorption cor- 
rected 2-10 keV luminosity of the broad-line AGNs are 
shown. In this figure, broad-line AGNs with R — i color 
redder than 0.3 are marked with open circles. Those 
with a R — i color of 0.3 arc redder than the scatter 
of typical broad- line AGNs in the re dshift range ob- 
served in SDSS (jRichards et all 120031 ). Hereafter, we 
designate broad-line AGNs with R — i redder (bluer) 
than 0.3 as red (blue) broad-line AGNs. In the dia- 
gram, blue broad-line AGNs follow the relationship ex- 
pected from the typical SEDs of broad-line AGNs as 
a function of intrinsic l uminosity which is shown with 
solid line in the figure (jMarconi et al.ll200l . The dot- 
ted line in the figure shows the relationship determined 
with the BCES bisector analysis for the blue broad-line 
AGNs. In lMarconi et~aTl ||200l . the SEDs around 3000 A 
are described with a power-law with a = —0.44 with 
L v cx v a and a dependence on optical-to-X-ray flux 
ratio, apXi on opt ical luminosity of broad-line AGNs 
(|Vignali et al.|[2003l ) is considered. The -B-band luminos- 
ity used in lMarconi et al.l ()2004f ) is converted to a 3000 A 
monochromatic luminosity assuming a typical SED of 
broad-line QSOs (1 549 x A4360L A 4360 = A3000L A 3ooo; 
iRichards et al.ll2006D . The consistency of the distribu- 
tion with this relation suggests that the blue broad-line 
AGNs have SEDs consistent with the optical-to-X-ray lu- 



minosity ratio, a ox , dependence on luminosity for typical 
non-absorbed broad-line AGNs. 

On the contrary, most of the red broad-line AGNs have 
a systematically fainter 3000A monochromatic luminos- 
ity than the blue broad-line AGNs at the same absorption 
corrected 2-10 keV luminosity. The fainter 3000A lumi- 
nosity suggests that the red broad-line AGN are affected 
by mild dust absorption although most of them show a 
strong Mg II broad-line. Additionally, some of the red 
broad-line AGNs are brighter in their 3000A luminosity. 
Most of them have the lowest absorption-corrected 2- 
10 keV luminosity and their red colors can be explained 
by contamination by a host galaxy component in the 
wavelength range. In both cases, the 3000A luminos- 
ity is not a good indicator of intrinsic luminosity, thus 
we use absorption-corrected X-ray luminosity instead of 
the 3000A monochromatic luminosity for the black hole 
mass estimation. We convert the absorption-corrected 
hard X-ray luminosity of the red broad-line AGNs to the 
intrinsic 3000A mon ochromatic lum i nosity using the re- 
lationship derived bv lMarconi et all (|2004[ ). Considering 
the scatter of blue broad-line AGNs around the relation- 
ship, we estimate the rms uncertainty of the intrinsic 
3000A monochromatic luminosity is 0.54 dex, which cor- 
responds to 0.27 dex in the Mbh uncertainty. We use the 
3000A monochromatic luminosity derived from the hard 
X-ray luminosity of 59 red broad-line AGNs out of the 
215 broad-line AGNs for the Mbh estimate. There are 
26 red broad-line AGNs among the 116 broad-line AGNs 
in the redshift range between 1.18 and 1.68 where the 
broad-line AGN BHMF and ERDF are derived below. 

In summary, 3000A monochromatic luminosities are 
derived in three ways, from the power-law component 
of the fitting of optical spectrum, the optical broad-band 
photometry for broad-line AGNs only with the NIR spec- 
trum, and hard X-ray luminosity for mildly-obscured or 
less-luminous broad-line AGNs. The resulting 3000A 
monochromatic luminosities are summarized in Table [3] 

4. BLACK HOLE MASS, BOLOMETRIC LUMINOSITY AND 
EDDINGTON RATIO 

The distribution of broad-line AGNs in the measured 
Mg II FWHM vs intrinsic 3000A monochromatic lu- 
minosity plane is shown in the left panel of Figure HI 
Large open squares indicate the FWHMs that are esti- 
mated with the broad Ha line and Equation (2). In the 
panel, the constant Mbh line derived from Equation (1) 
is shown with solid lines. From bottom to top, the lines 
correspond to Mbh of 10 6 , 10 7 , 10 8 , 10 9 , 10 10 , and 10 11 
M . The SXDS broad- line AGNs cover the M B h range 
between 10 7 ~ 10 10 M Q with a median of 3.2 x 10 8 M o . 

Although we define broad-line AGNs as having 
FWHMs above 1000 km s _1 , a few objects have FWHMs 
between 1000 and 2000 km s _1 . In the panel, we compare 
the distribution with tha t of broad-line AGN s from SDSS 
DR5 (filled gray circles; [SIien~et~all l2008ab . The 49,526 
SDSS broad-line AGNs, which arc selected with broad- 
line component whose FWHM is larger than 1200 km 
s _1 , are distributed between z — 0.3 and 2.4. Though 
there are far larger number of broad-line AGNs in the 
SDSS sample than in the SXDS sample, again only a neg- 
ligible fraction of broad-line AGNs have FWHMs smaller 
than 2000 km s _1 . A similar rapid decrease of broad- line 
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Figure 8. Mg II FWHM vs. 3000A monochromatic luminosities of SXDS broad-linc AGNs (filled circles). Broad-line AGNs whose Mg II 
FWHMs are estimated based on Ha FWHM are marked with large open squares. The solid lines indicate constant Mbh values for A/bh 
of 10 7 , 10 s , 10 9 , 10 10 , and 10 11 Mq from bottom to top. The dotted lines show constant AEdd for logAsdd of 1, 0, —1, —2, and —3 
from bottom to top. The thick dotted line indicates log A sdd of 0- The thic k solid line shows Mg II FWHM of 1000 km s" 1 . SDSS DR5 
broad-line AGNs that have a Mg II FWHM measurement fShen et al. 2008a) arc plotted as gray circles. Right) Same for broad-line AGNs 
at z = 1.18 - 1.68 only. 
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Figure 9. Black hole mass versus Eddington ratio of broad-line 
AGNs at redshifts between 1.18 and 1.68. Broad-line AGNs whose 
Mg II FWHM are estimated from Ha FWHM are marked with 
large open squares. The dotted line shows the relationship between 
log Mbh and Agdd f° r a broad-line AGN with L2-10 kcV = 10 43 
erg s , corresponding to the faintest object in the sample. The 
thick solid line indicate the constant Mg II F WHM of 1000 km s -1 . 
The distribution of the SDSS DR5 sample l|Shen et al.M2008al l in 
the same redshift range is shown with the contours. 



AGNs with FWHM smaller than 2000 km s" 1 is als o 
reported bv lHao et all (|2005t) and iStern fc Laori (|2012t ). 
They select broad-line AGNs in the local Universe from 
SDSS galaxy as well as quasar samples with broad Ha 
line above 1000 km s" 1 . The distribution of Ha FWHMs 
shows a rapid decrease in the FWHM range below a few 
1000 km s _1 . The physical origin of the cut off in the 
FWHM distribution is unknown. 

In the right panel of the figure, only SXDS and SDSS 
broad-linc AGNs at redshifts between 1.18 and 1.68 arc 



plotted. The luminosity limits of the surveys define 
the left-hand envelopes of the distributions. It can be 
seen that the SXDS broad-linc AGNs cover a 3000A 
monochromatic luminosity down to 10 43 5 erg s _1 in the 
redshift range. This luminosity limit is more than an or- 
der of magnitude fainter than SDSS sample in the same 
redshift range. 

We estimate the bolometric luminosity, Lboi from 
the 3000A monochromatic luminosity and using a 
bolometric correction factor of 5.8 for the monochro- 
matic luminosity from Richards et al.l (|2006| ). following 
iVestergaard fc Osmerl (|2009j). In Table H M BH , L hoh 
and AEdd are tabulated. The bolometric luminosities 
of the SXDS sample are distributed between Lboi = 
10 45 - 10 46 erg s" 1 with a median of 2.9 x 10 45 erg 
s _1 . The Eddington ratio of each object is calculated as 
AEdd = ^boi/^Edd- ^Edd is the Eddington-limited lumi- 
nosity given by 1.26 x 10 38 Mbh (erg s _1 ) with Mbh in 
units of Mq. The dotted lines in Figure [8] represent con- 
stant log AEdd values with 1, 0, —1, —2, —3, and —4 from 
bottom to top. log AEdd of is shown with thick dotted 
line. The SXDS broad-line AGNs range in log AEdd be- 
tween and —2. 

The uncertainties of A/bh f° r broad-linc AGNs with 
FWHM of Mg II and 3000A monochromatic luminosity 
of the spectroscopic data are typically 0.1 dex from the 
FWHM measurement and 0.04 dex from the variability 
between the epochs of spectroscopic and imaging obser- 
vations. As we described in Section [3~2l there can be sys- 
tematic uncertainty of 0.08 dex due to the Fe II template 
uncertainty. For broad-line AGNs with only Ha mea- 
surement, the scatter between the FWHMs of Mg II and 
Ha and the uncertainty due to the 3000A monochromatic 
luminosity add Mbh uncertainties of 0.22 and 0.05 dex, 
respectively. Furthermore, for mildly-obscured or less- 
luminous broad-line AGNs whose 3000A monochromatic 
luminosity is estimated with hard X-ray luminosity, the 
additional Mbh uncertainty of 0.27 dex comes from the 
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Table 4 

Black hole mass and Agdd of broad-line AGNs a 



ID 


z 


Iog[M BH /M ] 


l°g[-£<Bol 
/ [erg s" 1 ] ] 


log ^Edd 


0010 


1.225 


8.8 


45.65 


-1.22 


0018 


1.452 


8.9 


45.92 


-1.12 


0019 


1.447 


8.9 


46.18 


-0.87 


0023 


1.534 


8.9 


45.87 


-1.16 


0027 


2.067 


9.0 


46.61 


-0.52 


0034 


0.952 


8.5 


45.49 


-1.11 


0036 


0.884 


8.5 


46.00 


-0.64 


0037 


1.202 


8.4 


45.21 


-1.30 


0050 


1.411 


7.9 


45.56 


-0.43 


0056 


1.260 


9.0 


45.50 


-1.60 



a Table 4 is published in its entirety in the electronic 
edition of ApJ. A portion is shown here for guidance re- 
garding its form and content. 



luminosity conversion. Because Lboi is determined with 
^A3000j which is used for the Mbh estimation, the uncer- 
tainty of the AEdd is the same for the uncertainty of the 
M BH - 

The distribution of broad-line AGNs at redshifts be- 
tween 1.18 and 1.68 in the AEdd and Mbh plane is shown 
in Figure^ Because of the flux limit of the survey, no ob- 
ject in the lower left corner on the plane can be detected. 
It should be noted that the flux limit runs diagonally. 
The dotted line in the figure shows the constant hard X- 
ray luminosity line with log-L2-iokev(c r g s" 1 ) = 43. For 
10 8 M Q and 10 7 M Q , we can detect objects with log AEdd 
down to —1.5 and —0.5, respectively. The distribution is 
compared with those from SDSS (contour). Thanks to 
the deep detection limit of SXDS, we can select objects 
with lower mass as well as lower Eddington ratio in the 
same redshift range. In the figure, the relation between 
AEdd and Mbh with constant FWHM of 1000 km s _1 is 
shown with the thick solid line. The cutoff in the distri- 
bution of FWHM seen in Figure [3] appears as a deficit 
of objects in the upper left region with Mbh < 10 8 (M Q ) 
and log AEdd > —0.5, although the region is above the 
detection limit. 

5. BLACK HOLE MASS & EDDINGTON RATIO 
DISTRIBUTION FUNCTIONS 

5.1. Binned broad-line AGN BHMF and ERDF with 
Vniax Method 

First we derive the binned BHMF and ERDF for the 
broad-line AG Ns between 1.18 < z < 1.68 using the 
V max method (|Avni fc Bahcalllll980f ). Detailed numbers 
for the sample are shown in Table [5J In the calcula- 
tions of the binned BHMF and ERDF, we only consider 
broad-line AGNs with AEdd larger than 0.01, and remove 
2 broad-line AGNs below the limit. We also remove 2 
broad-line AGNs (SXDS0613 and SXDS0738) with nei- 
ther Mg II nor Ha FWHM measurements in the redshift 
range. We derive the binned broad-line AGN BHMF and 
ERDF for the soft- and hard-band samples separately. 

For the binned broad-line AGN BHMF, we divide the 
mass range 6.6 < log ( Mbh /M©) < 10.2 into 9 bins with 
bin width, A log Mbh, of 0.4 dex. The number density 
in a bin between (log Mbh — A log Mbh/2) ~ (log Mbh + 
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Figure 10. Binned active BHMFs for broad-line AGNs derived 
with soft- (filled circles) and hard- (open circles) band samples. 
Lines indicate corrected BHMF determined using the Maximum 
likelihood method with four different combinations of functional 
forms of BHMF and ERDF. Solid lines are for the doublc-powcr- 
law BHMF and dotted lines are for the Schechter BHMF. Thick 
lines are determined with the log-normal ERDF and thin lines are 
determined with the Schechter ERDF. Because the difference of 
the corrected BHMF with the different ERDFs is small, thick and 
thin lines overlap. 
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Figure 11. Binned ERDFs for broad- line AGNs derived with 
soft- (filled circles) and hard- (open circles) band samples. Lines 
indicate corrected ERDF determined using the Maximum likeli- 
hood method with four different combinations of functional forms 
of BHMF and ERDF. Solid lines are for the double-power-law 
BHMF and dotted lines are for the Schechter BHMF. Thick lines 
are for the log-normal ERDF and thin lines are for the Schechter 
ERDF. Because the difference of the corrected ERDFs with the 
different BHMFs is small, solid and dotted lines overlap. All of 
the ERDFs are normalized by matching the number density in the 
range logAEdd = —2.0 — 1.0 to the number density of broad-line 
AGN BHMF in the range logM BH = 7.0 - 11.0. 

A log Mbh/2) is given by, 

n 1 

< I ) bh(Mbh)A log Mbh = T7~- ( 3 ) 

i=i Va ' i 

The summation is done for the n broad-line AGNs in 
the mass bin. i is the index for a broad-line AGN in 
the mass bin. V a ,i is the effective survey volume for the 
i-th broad-line AGN in the comoving coordinate and its 
inverse represents the contribution of the broad- line AGN 
to the comoving number density of the mass bin. V a .i is 
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Table 5 

Binned broad-line AGN BHMF 



<log[A/ B H/M ]> a $ soft .(logA/BH) b N so{t * hard (logM B H) b iV hard 



7.2 


1.65 


+ 


3.80a 
1.36 


1 











7.6 


5.17 


± 


3.41 


1 


1.50 


± 


1.07 


2 


8.0 


11.60 


± 


2.43 


23 


27.46 


± 


15.50 


17 


8.4 


19.90 


± 


3.55 


36 


43.21 


± 


21.54 


30 


8.8 


18.62 


± 


3.36 


34 


19.93 


± 


5.21 


28 


9.2 


4.74 


± 


1.50 


10 


4.53 


± 


1.51 


9 


9.6 




















10.0 


0.47 


+ 


1.09c 
0.39 


1 


0.47 


+ 


1.09c 
0.39 


1 



a The central value of Mbh in each bin. The bin size is 0.4 dex and extends ± 

0.2dex from the central value. 

b In units of 10~ 6 ^Mpc" 3 (A log[A/BH/A^o]) _1 ] 

c The upper and lower limits are determined following Gehrels (1986). 



Table 6 

Binned ERDF at z = 1.43 



<logA Edd ) a $ soft (logA Edd ) b N solt *hard (logA Bdd ) b JVhard 



-2.00 


5.79 


± 


-1.75 


5.68 


± 


-1.50 


16.51 


± 


-1.25 


18.36 


± 


-1.00 


18.15 


± 


-0.75 


15.54 


± 


-0.50 


7.12 


± 


-0.25 


3.04 


± 


0.00 


0.76 


+ 


0.25 


0.76 


+ 



2.60 


6 


6.03 


2.43 


6 


6.12 


4.11 


18 


21.29 


4.53 


20 


20.80 


4.47 


20 


46.98 


3.77 


18 


15.23 


2.38 


10 


3.79 


1.52 


1 


4.95 


1.75c 
0.63 


1 


2.95 


1.74c 
0.63 


1 


0.76 



± 


3.49 


•1 


± 


2.89 


5 


± 


7.06 


11 


± 


7.77 


15 


± 


33.71 


16 


± 


4.16 


15 


± 


1.94 


5 


± 


2.28 


5 


± 


1.76 


3 


+ 


1.74c 


1 




0.63 



The central value of A Edd m each bin. 
In units of 10 -6 klpc" 3 ( A log[A Edd ])~ 



The upper and lower limits are determined following Gehrels (1986). 
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Figure 12. Binned ERDFs for soft-band broad-line AGNs in the 
two mass ranges lo 80-8 - 5 and if) 8 ' 5-90 Mq are shown with filled 
circles and filled diamonds, respectively. 



given by, 



V a ,i 



tt(L xl ,logN m . 



1 + z' 

1 + z CP . 



represent the rcdshift range for the binned 
broad-line AGN BHMF (1.18 and 1.68, respectively). 



f2(Lxi, log Nm, z 1 ) is the survey area which is calculated 
assuming that the i-th broad-line AGN observed at Zi 
with absorption-corrected luminosity Lxi, absorption hy- 
drogen column density log Nm is at rcdshift z'. With an 
estimated absorption-corrected Lxi and log Nm as de- 
scribed in Section [2j we calculate the predicted count- 
rate for each broad-line AGN with z' instead of the ob- 
served Zi. The same X-ray spectral model of AGNs is 
used, as explained in Section O The survey area for the 
sample with likelihood larger than 7 is determined as a 
function of count-rate for the o verlapping region of the 
X-ray and deep optical surveys (jUeda et al.l 120081 ). The 
logN-logS relation derived with the area curve for the 
SXDS X-ray sources is consistent wi th the relation de - 
tcrmined in deeper Chandra surveys (|Ueda et al.ll2008l) . 
The consistency implies the position-dependent detection 
limit of the X-ray survey is reproduced well in the esti- 
mated area curve. If the predicted count-rate of an object 
at a certain redshift z' is below the smallest count-rate 
limit of the survey, f2(Lxi, log Nm, z') becomes above 
that redshift. The factor ((14- z') / (1 + z ccn )) k corrects for 
the number density evolution with (l4-;z) fc within the red- 
shift range to determine the corresponding number den- 
sity at z cen of 1.43. However, the correction factor is neg- 
ligible even if we introduce rather strong number density 
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evolution with k = 4 — 5 observed for the X-ray luminos- 
ity fu nction of AGNs (jUeda et all 12003 lHasi'nger et~aTl 
I2005H because the redshift range for this calculation is 
narrow. Therefore, we neglect this term hereafter and fix 
the value k as for simplicity. Finally V a j is obtained 
by integrating the corresponding survey area for the ob- 
ject at redshift z' multiplied by the cosmological volume 
element (dV / dz')dz' of unit solid angle in the redshift 
range. The uncertainty of the binned broad-line AGN 
BHMF is estimated using Poisson statistics as, 



E 

,fe=i 



i 



V a ,i A log Mbh 



1/2 



(5) 



The resulting binned broad-line AGN BHMFs arc shown 
in Figure (TU] and Table [5] The filled and open cir- 
cles represent the binned BHMFs of soft- and hard- 
band broad-line AGN samples, respectively. Both of the 
binned broad-line AGN BHMFs peak at around Mbh of 
10 8 ' 5 M Q . The binned broad- line AGN BHMFs are con- 
sistent each other within the 1 a uncertainty. 

The binned broad-line AGN ERDF is derived in the 
same way for the binned BHMF using the V max method. 
We bin the sample in AEdd by dividing the AEdd range of 
—2.125 < log AEdd < 0.375 into 10 bins with a bin width 
AAEdd of 0.25 dcx. The number density of a certain 
Eddington ratio bin between ( log AEdd — A log AEdd/2) ~ 
(logA E dd + AlogA E dd/2) is given by, 



$A(A E dd)Alog A Ed d : 



n 

v as 



2 = 1 



(0) 



n is the number of broad-line AGNs in the Eddington 
ratio bin. V a ,i is the same effective volume for the i-th 
broad-line AGN used in the binned BHMF. The uncer- 
tainty for each Eddington ratio bin is again given by 



1 



^ VK.iAlogA 



Edd 



1 1/2 



(7) 



The resulting binned ERDF is shown in Figure [TT] and 
tabulated in Tabled The filled and open circles in the 
figure are the binned ERDFs for the soft- and hard- band 
samples, respectively. The binned ERDFs are consistent 
with each other within the 1 a uncertainty. In the binned 
ERDFs, broad-line AGNs in the entire mass range are 
considered. Due to the flux limit of the survey, the sam- 
ple covers only broad-line AGNs with larger Mbh in the 
lower AEdd bins. Therefore the shape of the binned ERDF 
can be affected by the flux limit. 

In order to examine the Mbh dependence of the binned 
ERDF, in Figure fl2l the binned ERDFs derived with the 
soft-band selected broad-line AGNs in the Mbh range 
between 10 8 -°- 8 - 5 M Q and that between 10 8 - 5 - 9 M Q are 
plotted with filled circles and filled diamonds, respec- 
tively. The overall shapes of the binned ERDFs do not 
significantly differ from each other within the uncer- 
tainty, except for the lowest log AEdd range where the 
lower mass binned ERDF can be affected by the flux 
limit. Due to the limited number and mass range of the 
SXDS sample, there is no signature of the dependence of 
the ERDF on M BH - 
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Figure 13. Effective volume of the SXDS soft-band sample in 
the redshift range between 1.18 and 1.68 as a function of 2-10 keV 
hard X-ray luminosity. 



5.2. Corrected broad-line AGN BHMF and ERDF with 
Maximum Likelihood Method 

Both the binned broad- line AGN BHMF and ERDF arc 
affected by the detection limit determined by the X-ray 
count rate; at the low-mass end of the binned BHMF the 
sample covers only high AEdd broad-line AGNs and at the 
low Eddington ratio end of the binned ERDF the sample 
does not include broad- line AGNs with low Mbh- Such 
detection limits are not corrected for in the calculations 
of the binned broad-line AGN BHMF and ERDF. 

The effects of the detection limit can be corrected 
through st atistical methods a ssuming the forms of both 
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lihood method to a sample of low-redshift broad-line 
AGNs detected in the ESO/Hamburg survey, assuming 
that the intrinsic ERDF does not depend on Mrh- On 
the other hand. IKellv et all ff2T)10h . IS hen fc Kellvl (|20T2h 
and iKellv fc Shenl (|2012f T apply a Bayesian approach 
(IKellv et all 12009ft to SDSS broad-line AGNs. They 
introduce a dependence of ERDF on the black hole mass. 
Furthermore, the statistical scatter in the virial black 
hole mass estimate is considered in the calculations; the 
scatter can broaden a peak in the intrinsic broad-line 
AGN BHMF and a steep slope of intrinsic BHMF at the 
high-mass end can become flatter in the binned BHMF. 
If the intrinsic BHMF has a peak at a certain black hole 
mass and there is a turn-over at the low-mass end, the 
scatter can also affect the binned BHMF in the low-mass 
end. 

Here, we apply Maximu m likelihood method used in 
iSchulze fc Wisotzkil (|20Tot) for the broad-line AGNs in 
the soft-band sample. Because there is no significant dif- 
ference between the binned broad-line AGN ERDFs in 
the high- and low-mass ranges as shown in Figure 1121 
we assume that the intrinsic broad-line AGN ERDF is 
constant regardless of black hole mass. The effect of the 
scatter of the virial black hole mass estimate is not con- 
sidered in this paper, because the SXDS sample covers 
smaller black hole mass and Eddington ratio ranges than 
the SDSS sample, most of the SXDS broad-line AGNs lie 
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Table 7 

Best fit parameters for the corrected broad-line AGN BHMF and ERDF 



BHMF a 


ERDF b 


(j>* c 


logM* d 


a 


/3 


Io g A Bdd e 




2DKS g Note h 


DPL 
SCH 
DPL 

SCH 


SCH 
SCH 
LOG 
LOG 


19.82j£» 
20.80iy° 

±».o»_ 3 20 
2f) fil + 3 - 30 


o H0.08 

8 -''-0.08 

8.59l° ii 
a 7B+0-08 
8 -' 8 -0.08 


U - UO -0.28 
u " 3o -0.14 
u - u9 -0 .27 
u -°°-0.14 


o ^7+0.28 
z -°'-0.39 

9 7n +0.29 
z -' u -0.40 


-0.581°,;" 
-0.59±g;iJ 

9 07+0. 16 
9 ao+0.16 


-0.74+g : » 

u -'^-0.15 
1 9^+0.14 

1 24+ ' 13 
• L-x -0.10 


86 

65 x 

79 

62 



a DPL: doublc-power-law, SCH: Schechter functions. 

b SCH: Schechter, LOG: log-normal distributions. 

c In unit of 10 _ 6 Mpc -3 per unit log Mbh bin. 

d In unit of M Q . 

e )x for log-normal distribution. 

f (7 for log-normal distribution. 

s 2DKS probability in unit of %. 

h Consistency with high-luminosity end of hard X-ray luminosity function. In details, sec Section 5.3. 



in the mass range below the knee of the BHMF, and in 
order not to be affected by the uncertainty associated 
with the modeling of the scatter. The high-mass end of 
the corrected broad-line AGN BHMF can be affected by 
the flattening due to the scatter. 

In this evaluation, we assume the shape of the corrected 
broad-linc AGN BHMF to be cither a doublc-power-law 
or a Schechter function, 



$bh(M bh ) 
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respectively, with the functions expressed per log Mbh- 
We introduce these two forms because the double-power- 
law describes the AGN luminosity function well, and 
the Schechter function describes the luminosity and mass 
functions of galaxies. Even though a modified Schechter 
functi on describes the non-active BHMF in the local Uni- 
verse ([Aller fc Richstond [20021: iShankar et al.ll2004D . we 
do not implement such a function, since it does not con- 
verge and results in a rather small M* (1O 5_6 M0). 

For the corrected ERDF we assume the Schechter func- 
tion and the log-normal distributions. 
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We do not consider the Mbh dependence of the ERDF. 
For both distributions, we define the normalized cor- 
rected ERDF with 



-P\ (AEdd) 



$(A E dd) 



$(A E dd)dlog A E dd 



(12) 



We normalize the corrected ERDF in the range log A E dd 
of -2.0 and 1.0. In these functions, <j>* , a, /J, and M* 
of the corrected BHMF and a\ and A Edd of the corrected 
ERDF (or a and fi for log-normal ERDF) are free param- 
eters. We derive the best fit parameters with th e Max- 
imum Likelihood method ((Marshall et al.l 11983ft . The 
likelihood function is written as 



N 



S-- 



-2y^ln p(M B H,i, A B dd,i) 
i=i 

+ 2 J /p(M B H,A E dd) dlogA Edd d log A/bh (13) 

and the model parameters that minim ize S are the best 
fit parameters ([Marshall et al.lll983D . The sum of the 
first term will be taken for the entire N broad-linc 
AGNs in the sample. The term j>(Mbh, AEdd) is the ex- 
pected number of black holes with Mbh and AEdd in 
unit log Mbh and log AEdd intervals in the survey red- 
shift range with the assumed 4>bh(Mbh) and Pa (AEdd) • 
p(M BH , A E dd) is given by 

p(Mbh, A E dd) 



^'(M B H,AEdd,z') 

Pa (AEdd) ^bh (Mbh) (^^-) ^z'. 

\l + z ccn J dz' 

(14) 

fi' (Mbh, AEdd j is the survey area for an broad-linc 
AGN with Mbh and A E dd at z' . We calculate the ex- 
pected count rate for each combination of Mbh, AEdd, 
and z' assuming the same model for the X-ray spectrum 
of AGNs used in Section [2] We convert the Lboi for 
a combination of Mbh and AEdd to L2-10 keV with th e 
bolometric correction factor from iMarconi et ail ([2004D . 
In this calculation we do not consider the effect of X-ray 
absorption because the effect is not significant for the 
sample of broad-line AGNs. In Figure [T3l we plot the 
effective volume 



cir 



n' (M BH , A Bd d, ( 15 ) 



as a function of -L2-10 kcV- We consider log AEdd down 
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to —2.0 and Mbh of 10 6 M Q . Again we use k of as 
explained in the previous subsection. We minimize S 
with the 6 free parameters with the downhill simplex al- 
gorithm (jNelder fc Meadlll965| ). The la uncertainty of 
the best-fit parameters for each model can be evaluated 
by the increase of S by 1 from the minimum value. In 
order to determine the la uncertainty of each parame- 
ter, the parameter is changed from its best value to a 
different value, and fixing the parameter at the value, 
the same minimization process is applied for the other 
parameters and we evaluate the change of the minimum 
S value from the best fit value of S. The uncertainty of 
the selected parameter is determined by the change of 
the minimum of 1 from the best fit value. The resulting 
best-fit parameters are summarized in Table [3 

The resulting corrected broad-line AGN BHMFs and 
ERDFs are shown in Figures [TU] and [TTJ The corrected 
ERDFs are normalized by matching the number density 
in the range logAEdd = —2.0 — 1.0 and that of the cor- 
rectedBRMF in the range log M B h = 7.0-11.0. The cor- 
rected BHMF follows the estimated number density with 
the Vmax method well in the mass range above 10 85 M Q . 
This is consistent with the limit of the SXDS sample; it 
extends down to a Eddington ratio of 0.01 in the mass 
range above ~ 10 8 Mq. In the lower mass range, the cor- 
rected BHMF is slightly larger than the number density 
derived with the VJ nax method. The estimated correction 
is consistent with the corrected ERDF; for 10 7 ' 5 M Q black 
holes the sample covers a Eddington ratio of ~ 0.1, and 
the ratio below and above this Eddington ratio is 1.7—1.8 
from the corrected ERDF. Therefore, the estimated cor- 
rection is not large. The shape of the corrected ERDF 
is consistent with that derived in the 10 s — 10 9 M Q mass 
range with the V max method. The correction required for 
the ERDF is only significant in the low Eddington range, 
logAEdd below —1.5. Thanks to the coverage for rela- 
tively low-luminosity broad-line AGNs, the SXDS sample 
is rather complete in a wide mass and Eddington ratio 
range. 

The minimum value of the likelihood function does not 
reflect the goodness of the fit. We evaluate the goodness 
of the fit by applying the two-dimensional Kolmogorov- 
Smirnov (2DKS) test on t he distribution of the sample i n 
the Mbh and AEdd plane (|Fasano fe Franceschini|[l987l ). 
The resulting 2DKS probability for the deviation in the 
plane are shown in Table [7] All of them exceed 20% and 
all of the models fit well the distribution of objects on 
the plane. 

In the mass range below 10 8 ' 5 M Q , the corrected 
BHMFs show possible decline to the low-mass end. The 
mass range is well above the detection limit as shown 
in Figure [5J but the decline can be affected by the com- 
pleteness of the broad-line AGN sample. For example, 
there is a possibility that among objects without spec- 
troscopic identification low-luminosity broad-line AGNs 
whose SEDs are dominated by host galaxy components 
are classified as narrow-line AGNs in the photometric 
rcdshift estimation, and they would be missed in the cur- 
rent broad-line AGN sample. 

5.3. Constraint from the Hard X-ray Luminosity 
Function 

The luminosity function of broad-line AGNs is the con- 
volution of their BHMF and ERDF, therefore we can con- 



strain the shapes of broad-line AGN BHMF and ERDF 
further by using the luminosity function determined from 
a combination of various AGN samples. In particular, 
the number density at the bright end of the luminos- 
ity function obtained through wider but shallower sur- 
veys than SXDS can constrain the shapes of broad-line 
AGN BHMF and ERDF in the high M BH and A Ed d 
range. It needs to be noted that the fraction of ob- 
scured narrow-l ine AGNs is low in the luminosity range 
(jHasin gor 20081) and the luminous end of the luminosity 
function is thought to be dominated by broad-line AGNs. 

In Figure 1141 we plot the observed hard X-ray lumi- 
nosity function of AGNs at z = 1.2 — 1.6 based on a 
combined sample of AGNs from various hard X-ray sur- 
veys as filled squares (Ueda et al. 2012, in preparation). 
The horizontal axis is the absorption-corrected 2-10 keV 
luminosity. The hard X-ray luminosity functions of the 
broad-line + narrow-line AGNs in the SXDS soft-band- 
selected and hard-band-selected samples are shown as 
open circles and squares, respectively. The hard-band- 
selected AGN luminosity function in the SXDS is consis- 
tent with that derived from the combined sample. The 
soft-band-selected AGN luminosity function has slightly 
lower number density than the hard-band-selected AGN 
luminosity function below £2-10 koV = lO^ergs -1 . The 
lower number density can be qualitatively explained with 
the fact that the soft-band-selected sample misses heav- 
ily obscured AGNs and the fraction of obscured AGNs 
is larger for AGNs with lower luminosity. Because all 
of the luminosity functions are consistent each other in 
the luminosity range above L2-10 kcV = 10 44 erg s _1 , 
therefore we use the number density of AGNs at the lu- 
minous end of the combined sample as the constraint on 
the broad-line AGN BHMF and ERDF. 

The filled circles represent the hard X-ray luminos- 
ity function of broad-line AGNs in the soft-band sam- 
ples. Again the hard X-ray luminosity is corrected for 
intrinsic absorption. The hard X-ray luminosity func- 
tion shows turn over at L2-10 kcV = 10 44 (erg s _1 ). 
Such turn over is also observed in the broad-line AGN 
hard X-ray l uminosity function a t z = 1 — 3 from Chan- 
dra surveys (jYencho et al.l l2009t ) . At least part of the 
decrease can be explained with the increasing fraction 
of obscured narrow-line AGNs in the lower luminosity 
range as described below. The lines in the figure show the 
hard X-ray luminosity function derived from the convolu- 
tion of the corrected broad-line BHMF and ERDF. Solid 
and dotted lines correspond to luminosity functions with 
the corrected BHMF of doublc-power-law and Schcchtcr 
forms and the thick and thin lines are those derived with 
Schechter and log-normal ERDF, respectively. The hard 
X-ray luminosity for each set of Mbh and AEdd is derived 
with the bolometric correction for 3000A monochromatic 
luminosity and the relation between A3000La3ooo and 
^2-iokcV shown with the solid line in Figure [7] Con- 
sidering the upper limit of the obser ved Mbh of broad- 
line AGNs (IVestergaard et al.1 120081 ) and local galaxies 
(jMcConnell et al.ll201ll ). we limit the M B h range to be 
10 6-io.5 Mq and the i og A Edd range to be -2 and 1. 

The hard X-ray luminosity functions derived with the 
double-power-law BHMF or log-normal ERDF can repro- 
duce the observed number density at the high-luminosity 
end. On the contrary, if both the BHMF and ERDF are 
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Figure 14. Left) Hard X-ray luminosity function of X-ray-selected AGNs (Ueda et al. 2012, in preparation; filled squares) compared 
with that reproduced from the corrected broad-line AGN BHMFs and ERDFs (lines). The open circles and squares show the luminosity 
function derived from the soft-band-selected and hard- band-selected AGN samples of SXDS in the range 1.18 < z < 1.68, respectively. 
Filled circles represent the luminosity function of broad-line AGNs in the soft-band sample, which is used to derive the corrected broad-line 
AGN BHMF and ERDF. The hard X-ray luminosities of all AGNs are corrected for intrinsic absorption. Solid and dotted lines correspond 
to the luminosity functions for the BHMF characterized by doublc-powcr-law and Schechter functions and the thick and thin lines are 
those derived with the Schechter and log-normal ERDF, respectively. The thin dot-dashed line shows the luminosity function derived with 
double-power-law BHMF + log-normal ERDF with limiting Mbh range below 10 10 Mq. The thin long-dashed line represents the luminosity 
function with doublc-powcr-law BHMF + log- normal ERDF after correcting for the obscured fraction following Hasingcr (2008). Right) 
Zoom in version of the left panel in the high-luminosity range. 



modeled with an exponential-cutoff such as the Schechter 
BHMF with Schechter ERDF (thin dotted line), the high- 
luminosity end of the luminosity function cannot be re- 
produced; the predicted number density is more than 
one order of magnitude smaller than the observed lu- 
minosity function. Therefore such models are unlikely 
to represent the BHMF and ERDF of broad-line AGNs 
at z = 1.4. Furthermore, if we limit the mass range of 
BHMF up to 1O 1O M , the predicted luminosity function 
with a double-power-law BHMF and log-normal ERDF 
are represented by the dot-dashed line in the figure; the 
predicted density in the high luminosity end is much 
lower than the observed one. Thus in order to repro- 
duce the number density of luminous AGNs, the BHMF 
needs to be extended up to 10 10 ' 5 Af Q under the assump- 
tion that the ERDF is constant over the wide Mbh range. 

In the luminosity range below L2_io koV = 10 44 erg 
s _1 , the broad-line AGN luminosity function shows de- 
cline toward lower luminosity. At least part of the de- 
cline can be explained with the increasing fraction of ob- 
scured narrow-line AGNs in the lower luminosity range. 
In Figure [T4l we also plot model luminosity function de- 
rived with a double-power-law BHMF and log-normal 
ERDF af ter correc t ing th e fraction of narrow-line AGN 
following lHasingerl (|2008l ) with thin long dashed line. In 
IHasingeii (|2008D . luminosity dependence of the fraction 
is derived as a function of redshift. We use the rela- 
tion derived at z = 1.2 — 1.6. Although the corrected 
luminosity function still has lower number density than 
the broad-line + narrow-line AGN luminosity function, 
the large fraction of the discrepancy between the broad- 
line luminosity function and total luminosity function 



seems to be explained with the obscured fraction down 

to £2-10 kcV = 

10 43 erg s _1 . The remaining discrep- 
ancy can be caused by the possible incompleteness of 
the broad-line AGN BHMF due to the spectroscopic in- 
completeness and the difficulties in identifying broad-line 
AGNs with substantial host contamination. The uncer- 
tainty of the fraction of narrow-line AGN is still large and 
the sample size of the SXDS is limited. Larger sample is 
necessary to understand the remaining discrepancy. 

6. DISCUSSION 

6.1. Broad-line AGN BHMF at z and Its evolution 

to z = 

The binned and corrected broad-line AGN BHMFs 
from the soft-band sample are compared with the binned 
and esti mated broad-l i ne A GN BHMF from the SDSS 
sample (|Shen fc Kellvl [20T1 in Figure Q21 The SXDS 
binned BHMFs (filled circles) is larger than the binned 
BHMF of the SDSS sample (open circles) by 1 order of 
magnitude at 1O 8 ' 5 M , and by two orders of magnitude 
at 10 8 Mq. This is mostly because broad-line AGNs with 
lower Eddington ratio are detected in the SXDS sample 
than in the SDSS; the SDSS sample is only 30% com- 
plete do wn to Mbh of ~ 1 9 M^ and AEdd of ~ 0.6 at 
z = 1.4 (Kell y et all 120101 ). Additionally, as discussed 
in Section 13.41 the SXDS sample covers even mildly ob- 
scured broad-line AGNs that may be missed in the SDSS 
selection of broad-line AGNs due to color and stellarity 
issues. The thin solid lines in the figure represent the up- 
per and lower envelopes of the esti mated broad-line AGN 
BHMF with a Bayesian approach (jShen fc Kellvl 12012ft . 
The estimated number density at 10 8 M Q is consistent 
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Figure 15. Binned (filled circles) and corrected broad-line AGN BHMFs (thick solid line with double-power-law model and thick dashed 
line with Schechter model) from SXDS compared with binned (open circles) and estimated (thin solid lines) broad-line AGN BHMFs at 
z ~ 1.4 from SDSS IjShen & Kcllv 20l3). Upper and lower thin solid lines enclo s e the 68% possible area. The thin dashed line shows the 
corrected b road-line AGN BHM F in the local Universe from Schulzc & Wisotzki (201JJ). The thin dotted line indicates the total BHMF at 
z = 6 from lWillott et al.1 ( |20Tot ). 



with the binned and corrected BHMF of SXDS. 

We examine the evolution of broad-line AGN BHMF 
from z = 1.4 to by comparing the z ~ 1.4 cor- 
rected BHMF with that in the local Universe. The 
dashed lines in the figure show the corrected local ( z < 
0.3 ) broad-line AGN BHMF from iSchulze fc Wisotzkil 
(|2010T) . The binned BH MFs of broa d- line AGNs 
in the local Univer s e of iGreene fc Hoi (|2009l ) and 
iVestergaard fe Osmeri (120091) are co n sisten t with the 
binned BEMF of ISchulze fe Wisotzkil (f2010h . The cor- 
rected BHMF at z ~ 1.4 exceeds that of local Universe 
in the mass range above 10 s M@. However the behav- 
ior in the lower mass r ange is different; the local broad- 
line AGN BHMF from ISchulze fe Wisotzkil pOlOh shows 
a steep increase with decreasing black hole mass even 
below 10 7 M Q mass range. On the contrary, the SXDS 
sample shows a hint of turn-over at a mass of 1O 8 ' 5 M0. 
The difference in the corrected broad-line AGN BHMF 
may be indicative of a down-sizing trend of accretion ac- 
tivity among the SMBH population. As mentioned in 
Section 15.21 it needs to be noted that the identification 
of broad-line AGNs could be incomplete in the low-mass 
end. 

It is also possible that decline of the activity in the low- 
mass range from z = to 1.4 is caused by luminosity and 
rcdshift dependent obscuration to the nucleus. Although 
broad-line AGNs have a wide range of AEdd, the fraction 
of broad-line AGNs with low-mass SMBH is higher in 
the lower luminosity range. A luminosity dependence of 
the obsc ured fraction is observed in various AGN sam- 
ples (e.g. lAkivama et al.ll2000HUeda et aO2003t iSimpsonl 



120051 : iBrusa et al.ll2010l ); the obscured fraction is higher 
among lower luminosity AGNs. Recently, based on large 
samples of X-ray selected AGNs, it has also been sug- 
gested that the fraction of obscured AGNs is higher a t 
higher redshifts (e.g. iHasingerl l200l Iffiroi et al.N2012[ ). 
Such redshift and luminosity dependent obscuration may 
hide activity among low-mass SMBHs at high-redshifts. 
We examine the contribution from obscured narrow-line 
AGNs to the total active BHMF in Section IQ1 

6.2. Comparison with total BHMF at z ~6 

In order to examine the growth of SMBHs at higher 
redshift, we compare binned and correct ed broad-line 
AGN BHMF at z = 1.4 with that at z = 6. iWillott et all 
(|2010h examined the total BHMF at z ~ 6 using the op- 
tical luminosity function of broad-line AGNs at that red- 
shift. They derive the total BHMF from the luminosity 
function using an Eddington ratio distribution derived 
from a part of their sample. They estimate the total 
(active + non-active) BHMF by crudely correcting the 
obscured fraction and duty cycle (active fraction). The 
resulting best-estimate total BHMF is shown with the 
thin dotted line in FigureHS] They argue the total BHMF 
is constrained down to 10 8 Af Q . At that mass, the z = 6 
BHMF has a density of 2.5 x 10~ 7 Mpc~ 3 . The number 
density matches the number density of z = 1.4 corrected 
BHMF of broad-line AGNs at 10 9 5 M Q . If we naively 
assume that the 10 8 M Q black holes at z = 6 grow to the 
10 9 5 Mq black holes at z = 1.4 through accretion, then 
this implies mass growth by a factor 30 between z = 6 
and 1.4, i.e., in 3.5 Gyr period. Such growth implies 
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that the multiple of the Eddington ratio and duty cycle 
has order of 0.1 in that period for the massive SMBHs, 
if we use the rest mass energy to radiation energy con- 
version efficiency of 0.1. It needs to be noted that the 
high-mass end of the corrected BHMF can be affected 
by the flattening due to the uncertainties of the virial 
black hole mass estimates, therefore the effect needs to 
be corrected before further quantitative evaluation of the 
growth of SMBHs between z = 6 to 1.4. 

6.3. Eddington Ratio Distribution Functions 

Both of the binned and corrected broad-line AGN 
ERDF show a decline at Eddington ratio of 1. Such a 
distribution suggests that the accretion of X-ray-selected 
broad-line AGNs is limited by the Eddington luminos- 
ity. In Figure 1161 we compare the corrected broad-line 
AGN ERDF with the estimated broad-line AGN ERDF 
at z ~ 1.4 from the SDSS sample [thin solid lines, Shen 
& Kelly 2012). The shapes of the ERDFs match rather 
well. It needs to be noted that there is a possibility 
that we miss broad-line AGNs accreting at a rate close 
to the Eddington limit because our X-ray selection is 
in relatively high energy range; the energy range of the 
soft-band sample corresponds to 1.2-4.8 keV at z = 1.4. 
Optical to X-ray SED models of broad-line AGNs predict 
weak X-ray emission with steep X-ray spectrum among 
broad -line AGNs accreting with close to the E ddington 
limit (|Kawaguchi et al.ll200H iDone et al.ll2012h . 

In the same figure, we also plot the corrected 
broad-line AGN ERDF in the local Universe from 
ISchulze fc Wisotzkil (|2010[) . The shape of their binned 
ERDF in the local universe is co n sisten t with that de- 
rived by iKauffma nn fc Heckmanl ([2009) for a narrow- 
line AGN sample from SDSS. The shape of the corrected 
ERDF in the local Universe has a knee at a similar Ed- 
dington ratio (log AEdd of —0.6; Schulze & Wisotzki 2010) 
to that of the corrected ERDF at z = 1.4, but shows a 
steeper increase in the low Eddington ratio range down 
to log AEdd of —2.0. The evolution of the corrected broad- 
line AGN ERDF from z= 1.4 to indicates the fraction 
of broad-line AGNs with Asdd close to 1 is higher at 
higher rcdshifts. 

6.4. Contribution from Obscured Narrow-line AGNs 

Among the X-ray-selected AGNs in the redshift range, 
more than half of the AGNs are obscured narrow-line 
AGNs as shown in TableO In this subsection, the contri- 
bution of these obscured narrow-line AGNs to the binned 
active BHMF is evaluated using the hard-band sample 
which is less biased against obscured AGNs than the 
soft-band sample. In the calculation of the binned ac- 
tive BHMF, not only spectroscopically- identified narrow- 
line AGNs but also narrow-line AGNs with a photomet- 
ric redshift are included. Because we use the sample 
of the hard X-ray selected AGNs, the contribution from 
the heavily-obscured Compton-thick AGNs are not in- 
cluded. The black hole mass of the obscured narrow-line 
AGNs cannot be estimated from the FWHM of broad- 
line. Here, their black hole mass is estimated assuming 
they have same AEdd as broad-line AGNs with the same 
bolometric luminosity. 

First, the relation between the Eddington ratio and 
the bolometric luminosity is determined for the broad- 
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Binned active BHMF including Obscured 
Narrow-line AGNs 
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a The central value of AfeH m each bin. The bin 
size is 0.4 dex and extends ± 0.2dcx from the cen- 
tral value. 

b In unit of 10" 6 [Mpc- 3 (A logfMBH/Ai©]) -1 ] 

c The upper and lower limits are determined fol- 
lowing Gchrels (1986). 

line AGNs in the SXDS sample. The Lboi vs. AEdd dis- 
tribution of the broad-line AGNs is shown in Figure [171 
There is a correlation between the Eddington ratio and 
the bolometric luminosity. The solid line in the figure 
represents the relation determined with the least square 
fitting, 



log A Edd = 0.469 x logL bol - 22.46. 



(16) 



The relation is mainly driven by the virial black hole mass 
estimator used, and the scatter reflects the distribution 
of the FWHM. The scatter of the AEdd from the best fit 
relation is 0.4 dex, which is roughly consistent with the 
width of the distribution of the FWHM. In the figure, 
the distribution of SDSS dShen et al.ll2008aD and Large 
Brigh t Quasar Survey (LBQS) (jVestergaard fc Osmerl 
120091 ) broad-line AGNs are also shown with the contours 
and crosses. For the SDSS sample, only broad-line AGNs 
with Mbh from Mg II measurements are shown. All of 
the broad-line AGNs follow the same trend of increas- 
ing lo g AEdd with increasing -Lb i as discussed in iCrooml 

poll . 

For obscured narrow- line AGNs, their bolometric lu- 
minosities arc estimated from the absorption corrected 
2-10 keV luminosity with a bolometric correctio n factor 
for the 2-10 keV luminosity (jMarconi et al.ll2004l ). We do 
not use UV or optical luminosities to estimate bolometric 
luminosity, because they are severely suffered from dust 
extinction and contamination from host galaxies. Using 
the above Eddington ratios for three bolometric luminos- 
ity bins, we estimate black hole mass of each narrow-line 
AGNs as 



M BH [M G 



i bo i[erg s" 



1.26 x 10 38 x A 



Edd, median 



(17) 



The scatter of the iboi and AEdd relation is large, thus 
the black hole mass estimation is only valid in statistical 
sense. The binned BHMF derived including the narrow- 
line AGNs is shown in Figure [15] and is tabulated in 
Table [Sj The contribution of the obscured narrow-line 
AGNs is large in the mass range, log Mbh < 8.5. 

In order to examine the fraction of active black holes in 
the entire SMBH population, we compare the binned ac- 
tive BHMF including narrow-line AGNs with non-active 
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Figure 16. Binned (filled circles) and corrected broad-linc AGN ERDF (thick solid line with Schccheter E RDF and thick dott ed line with 
log-normal ERDF) from SXDS compared with the estimated broad-line AGN ERDF at z ~ 1.4 from SDSS l lShen fc K elly 2012) (thin solid 
lines). The upper and lower thin solid lines represent the 68 % enclosed region. The thin dashed line shows the corrected broad- line AGN 
ERDF in the local Universe from lSchulze fc WisotzEl J20TOT) . 
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Figure 17. Bolometric luminosity versus Eddington ratio of 
SXDS broad-linc AGNs (filled circles). The solid line represents 
the least square fitting result for the SXDS AGNs. The distribu- 
tion of SDSS DR5 broad- line AGNs whose Mbh is estimated from 
Mg II |Shcn ct al. 20083) is shown b y the contours. Crosses rep - 
resent broad- line AGNs from LBQS {Vcstcrgaard & Osmcr1 l2009l ) . 



BHMF at intermediate rcdshifts derived from the galaxy 
luminosity functions and Mbh vs . -kbuige relation ship 
(jTamura et al.|[200l iLi et all [2^11 . iLi et all poTl es- 
timate the non-active BHMF up to z ~ 2 from the galaxy 
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Figure 18. Binned active BHMF with (open triangles) and with- 
out (open circles) obscured narrow-line AGNs derived from the 
hard-band sample. Solid, dotted, and dashed lines indicate the 
non-active BHMFs at z = 0, 1, and 2 derived with the it-band 
galaxy luminosity functions from ILi et all 1 120111) . 

if-band luminosity function, stellar mass function, and 
their redshift evolution. They assume an evolution of 
the Mbh vs. ibuige relationship as a function of red- 
shift of the form MB H/^snhnmid.K oc (1 + z) 1A following 
iBennert et all (|2010( ). They also consider the redshift 
evolution of the average bulgc-to-total luminosity ratio, 
even though they do not include luminosity dependence 
of the ratio. The resulting non-active BHM F is consis- 
tent w ith that estimated in the local Universe ()Vika et all 
l200l . In then gure, the estimated non-active BHMFs at 
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z = 0, 1, and 2 from lLi et al.l ([201 1[ ) are plotted with the 
solid line, dashed line, and dotted line, respectively. The 
binned active BHMF including narrow-line AGNs lies be- 
tween the non-active BHMFs at z = 1 and 2 in the mass 
range above Mbh = 1O 8 ' 5 M ; a rather high fraction of 
active SMBHs in the mass range is implied at z = 1.4. 
It needs to be noted that the binned active BHMF is 
not corrected for the scatter of the virial black hole mass 
estimate and that of the relation between the log Lboi 
and the log AEdd- Therefore, the number density can be 
overestimated in the high-mass end due to the contami- 
nation from the objects in the lower mass range. In order 
to discuss the fraction of active SMBHs quantitatively, 
the scatters needs to be corrected. 

7. SUMMARY 

Black hole masses of X-ray-selected broad-line AGNs 
detected in the SXDS are estimated from the width of 
the Mg II broad-line and the 3000A monochromatic lumi- 
nosity. Because optical spectroscopic observations cov- 
ering the Mg II wavelength range are not complete for 
the entire X-ray-selected broad-line AGNs, the width of 
Ha broad-line measured with NIR spectroscopic survey 
is also used to provide a supplementary estimate if black 
hole masses for broad-line AGNs whose Mg II broad- 
line width is not available. The sample of broad-line 
AGNs is selected using X-ray detection as well as de- 
tection of a broad emission line of either Mg II in the 
optical or Ha in the NIR. Some of the broad-line AGNs 
have red rest-frame UV-colors, suggesting mild obscu- 
ration to their nucleus. For such broad-line AGNs, the 
hard X-ray luminosity is used to estimate their intrinsic 
3000A monochromatic luminosity assuming typical opti- 
cal to X-ray luminosity ratio as a function of bolometric 
luminosity. In total, black hole masses are estimated for 
215 broad-line AGNs at redshifts between 0.5 and 2.3. 
In the redshift range between 1.18 and 1.68, the black 
hole mass estimate is highly complete thanks to the sup- 
plemental measurements of the Ha line width from the 
NIR spectra. 

Using the black hole masses of broad-line AGNs in 
the redshift range, binned broad-line AGN BHMF and 
ERDF are initially calculated using the V m&x method. 
The formed BHMF shows a peak at 1O 8 5 M . The binned 
ERDF has a steep decline at AEdd of 1 and a rather 
flat distribution below the Eddington limit. Because the 
sample is X-ray-flux limited, both the binned BHMF and 
ERDF are affected by the detection limit; the low Mbh 
end of the binned BHMF misses the low AEdd objects 
and the low AEdd end of the binned ERDF does not in- 
clude the low Mbh broad-line AGNs. The effect of the 
flux limit is corrected by assuming that the ERDF is 
constant regardless of the black hole mass. Applying 
the Maximum Likelihood method with appropriate func- 
tional shapes for the broad-line AGN BHMF and ERDF, 
we determine the corrected BHMF and ERDF. We do 
not correct for the effect of the uncertainties in the virial 
black hole mass estimates. Thanks to the faint detection 
limit, the sample extends down to AEdd of 0.01 in the red- 
shift range, the correction is rather small and the shapes 
of the corrected BHMF and ERDF do not significantly 
differ from those of binned BHMF and ERDF. 

The corrected broad-linc AGN BHMF peaks at 
1O 85 M and shows a possible decline at the low-mass 



end. The number density around 1O 8 M is consistent 
with that estimated from the SDSS sample. The shape 
of the corrected BHMF is completely different to that 
in the local Universe which shows a steep increase down 
to 1O 7 M . The evolution of the shape of the corrected 
BHMF of broad- line AGNs from z = 1.4 to z = may 
be indicative of a down-sizing trend of accretion activity 
among the SMBH population. 

The corrected broad-line AGN ERDF also shows a 
strong decline at the Eddington-limit and rather flat dis- 
tribution below the limit. The shape is also consistent 
with that estimated from the SDSS sample. The strong 
decline suggests that the accretion of X-ray-selected 
broad-line AGNs is limited by the Eddington luminos- 
ity. The shape of the corrected ERDF in the local Uni- 
verse has a knee at a similar Eddington ratio but shows 
a steeper increase in the low Eddington ratio range down 
to A Ed d of 0.01. The evolution of the ERDF from z = 1.4 
to z = indicates the fraction of broad-line AGNs with 
AEdd close to 1 is higher at higher redshifts. 

Using the hard-band sample of X-ray-selected AGNs, 
we evaluate the contribution of obscured narrow-line 
AGNs to the binned active BHMF. We estimate the black 
hole masses of narrow-line AGNs assuming the same Ed- 
dington ratio as a function of bolometric luminosity as for 
broad-line AGNs. Once the contribution from narrow- 
line AGNs is included, the binned active BHMF shows a 
comparable number density to the non-active BHMFs at 
z = 1 — 2 above Mbh of 1O 85 M . Such a high density of 
active SMBHs at z = 1.4 suggests that massive SMBHs 
have high active fraction in the redshift range. 
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